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PRÉFACE 



Les révolutions techniques les plus importantes et les plus riches de conséquences ne 
sont pas toujours celles qui sont les plus visibles pour l’utilisateur final du produit. 
Les méthodes modernes de traitement numérique du signal entrent dans la catégorie 
des révolutions techniques aux conséquences encore insuffisamment perçues et qui 
ne font pas la première page des journaux. 

Il est intéressant d’ailleurs de réfléchir quelques instants à la manière dont de 
telles techniques voient le jour. Le traitement par le calcul numérique d’un signal au 
sens le plus large du terme n’est certes pas en soi une idée nouvelle. Lorsque Kepler 
tirait les lois du mouvement des planètes des séries d’observations de son beau-père 
Tycho Brahé, c’est à un véritable traitement numérique du signal qu’il se livrait, le 
signal en l’occurrence étant constitué par les séries temporelles des observations de 
positions de Tycho Brahé. Mais ce n ’ est que dans le courant de ces toutes dernières 
décennies que le traitement numérique du signal est devenu une discipline en soi : 
c’est que la nouveauté tient à ce que l’on peut maintenant procéder, en temps réel, au 
traitement de signaux électriques, et ceci par des méthodes numériques. 

Pour que cette évolution soit possible, il fallait que des progrès techniques, dans 
de nombreux domaines, voient progressivement le jour, et tout d’abord, bien sûr, la 
possibilité d’acquérir, sous forme de signal électrique, des informations à traiter. 
Cela impliquait le développement progressif de tout ce qu’il est parfois convenu 
d’appeler les capteurs d’informations, lesquels peuvent aller, dans leur complexité, 
de la simple jauge de contrainte (mais il a fallu de nombreuses et difficiles recherches 
de physique des solides pour la rendre possible) au radar. 

Il fallait aussi que se développent, avec les prodigieux progrès de la micro- 
électronique, les outils technologiques permettant de réaliser, aux cadences extrême- 
ment élevées qu’implique le traitement en temps réel, des opérations arithmétiques 
que les premiers ordinateurs (l’ENIAC n’a que 40 ans) ne pouvaient réaliser qu’en 
plusieurs heures souvent interrompues de plusieurs pannes, et que nous trouvons 
aujourd’hui tout à fait naturel de voir exécutées par un micro-processeur de quelques 
grammes consommant seulement quelques milliwatts, et dont le temps moyen entre 
pannes dépasse dix ans. 

Il fallait enfin que les méthodes de programmation, c’est-à-dire d’utilisation 
optimisée de ces outils nouveaux, aient pu progresser, car quelles que soient les 
immenses capacités de calcul des micro-processeurs modernes, il n’est pas indifférent 
de ne pas gaspiller ces possibilités en opérations inutiles. L’invention des algorithmes 
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de transformée de Fourier rapide est un des exemples les plus frappants de cette 
importance des méthodes de programmation. Cette convergence des progrès tech- 
niques dans des domaines aussi différents relevant pour les uns de la physique, pour 
beaucoup de l’électronique, pour d’autres des mathématiques, n’a pas été acciden- 
telle. Dans une certaine mesure, chacun des progrès a suscité le besoin nouveau 
auquel un nouveau progrès dans un autre domaine permettait de répondre. Il serait 
sans doute utile, du point de vue de l’histoire et de l’épistémologie des Sciences et des 
Techniques d’entreprendre un jour une étude approfondie de ce cas. 

Car les conséquences en sont d’ores et déjà considérables. Sans doute le traite- 
ment analogique de signaux électriques a-t-il précédé le traitement numérique et sans 
doute continuera-t-il à occuper une place importante dans certaines applications, 
mais les avantages du traitement numérique qui tiennent en deux mots « précision et 
fiabilité » ont seuls rendu possibles certaines réalisations et qui déborden t de loin les 
secteurs de l’électronique et des télécommunications dans lesquels ces techniques ont 
vu le jour. Pour n’en citer qu’une, la tomodensitographie par rayon X, les «scanners » 
sont basés sur l’application d’un théorème dû à Radon et connu depuis 1917. Seules 
les évolutions que nous avons mentionnées plus haut ont rendu possible la réalisa- 
tion pratique de ce nouvel outil de diagnostic médical. Il y a gros à parier que les 
techniques de traitement numérique du signal trouveront demain leur place dans des 
produits de plus en plus variés, y compris les produits utilisés par le grand public 
qui, tout en bénéficiant des avantages de prix, de performance et de fiabilité que ces 
techniques rendent possibles, ne se rendra pas toujours compte de la prodigieuse 
imbrication de recherche, de technique et d’invention que suppose ce progrès. Cette 
évolution a d'ailleurs déjà commencé dans le cas des récepteurs de télévision. 

Mais lorsque se produisent de telles révolutions techniques, une autre difficulté 
se rencontre presque inévitablement. C’est celle de la formation des utilisateurs à ce 
qui est non seulement un nouvel outil, mais souvent un nouveau mode de pensée. 
Cette étape de la formation peut devenir, si l’on n’y prend garde, un goulot d’étran- 
glement dans l’introduction de nouvelles techniques. C’est pourquoi l’ouvrage de 
M. Bellanger, dont le point de départ est un enseignement donné depuis plusieurs 
années à l’École Nationale Supérieure des Télécommunications et à l’Institut 
Supérieur d’ Électronique de Paris, constitue un événement dont il convient de souli- 
gner l’importance. Ouvrage didactique, accompagné d’exercices, contenant plusieurs 
programmes, que certains pourront souvent utiliser tel quel, il contribuera sans 
aucun doute à accélérer encore une évolution désirable et nécessaire. 



P. AIGRAIN 
1981 
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AVANT-PROPOS 



L'innovation impose à l'ingénieur une mise à jour permanente de ses connais- 
sances et une bonne information sur le potentiel offert par les techniques nou- 
velles, découvertes et mises au point dans les laboratoires de recherche. En traite- 
ment du signal, les techniques numériques apportent des possibilités prodigieuses : 
la conception rigoureuse des systèmes, une grande reproductibilité des équipe- 
ments, une grande stabilité de leurs caractéristiques en exploitation et une remar- 
quable facilité de supervision. Cependant, ces techniques présentent un certain 
degré d'abstraction et leur application aux cas concrets requiert un ensemble de 
connaissances théoriques, jugées souvent plus familières ou plus facilement acces- 
sibles au chercheur qu'à l'ingénieur, et qui peuvent représenter un obstacle à leur 
utilisation. L'ambition du présent ouvrage est de vaincre cet obstacle et de faciliter 
l'accès aux techniques numériques en faisant la liaison entre la théorie et la pra- 
tique, et en mettant à la portée de l'ingénieur les résultats les plus utiles dans ce 
domaine. 

La base de cet ouvrage est un enseignement donné dans des écoles d'ingé- 
nieurs, d'abord l'École Nationale Supérieure des Télécommunications et l'Institut 
Supérieur d'Électronique de Paris, puis Supélec et le CNAM. Il s'adresse donc 
d'abord aux ingénieurs. L'auteur s’est efforcé d'y faire une présentation claire et 
concise des principales techniques de traitement numérique, de comparer leurs 
mérites et de donner les résultats les plus utiles sous une forme directement exploi- 
table aussi bien pour la conception des systèmes que pour une évaluation rapide 
dans le cadre de l'élaboration d'un projet en temps limité. Les développements 
théoriques ont été réduits à ce qui est nécessaire pour une bonne compréhension 
et une application correcte des résultats. Le lecteur trouvera dans les références 
bibliographiques les compléments qu'il pourrait souhaiter. A la fin de chaque cha- 
pitre, quelques exercices, souvent tirés de cas concrets, permettent de tester l'assi- 
milation de la matière du chapitre et de se familiariser avec son utilisation. Pour 
ces exercices, des éléments de réponse et des indications ont été regroupés en fin 
d'ouvrage. Il convient également de signaler que des efforts ont été faits pour 
introduire une terminologie française, qu’il serait souhaitable de compléter et 
généraliser afin de donner à notre langue sa place à part entière dans le domaine. 

Cet ouvrage s’adresse également aux chercheurs à qui il peut apporter, en plus 
d'un ensemble de résultats utiles, des indications pour l'orientation de leurs tra- 
vaux, en faisant clairement apparaître les contraintes de la réalité technique. Il 
contient de plus un certain nombre de résultats provenant des travaux de 
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recherche de l'auteur et de ses collaborateurs. En effet, pour établir le dialogue 
avec les chercheurs et être en état de faire bénéficier la technique de leurs décou- 
vertes dans les délais les plus brefs, l'ingénieur doit s'intégrer à la communauté 
scientifique et apporter sa propre contribution à la recherche ; par ses contacts per- 
manents avec les aspects concrets, il peut non seulement évaluer et conforter les 
résultats obtenus par les chercheurs mais encore ouvrir de nouvelles voies. 

Par rapport aux précédentes, cette huitième édition apporte des compléments, 
des parties nouvelles et des simplifications. Les compléments portent sur des résul- 
tats partiels nouveaux introduits dans différents chapitres, des expressions simpli- 
fiées et des éclaircissements. Les parties nouvelles se situent principalement aux 
chapitres 11 et 13. En effet, il est apparu intéressant de montrer les connexions 
entre les filtres QML et les ondelettes et que les ondelettes élémentaires résultent 
simplement d'un calcul de filtre avec des contraintes particulières. Quant à l'ana- 
lyse et modélisation des signaux, ces fonctions s'intégrent de plus en plus dans les 
systèmes. Le traitement du signal se généralisant avec l'électronique, le domaine 
d'application est de plus en plus vaste et il a semblé judicieux de se limiter à 
quelques illustrations des méthodes de base. 

Il faut souligner que les travaux sur lesquels est basé le présent ouvrage ont 
été à l'origine menés en collaboration et avec le soutien du Centre National 
d'Études des Télécommunications, à qui l'auteur tient à exprimer sa reconnais- 
sance. Il tient également à exprimer sa profonde gratitude à Monsieur J. Daguet, 
Directeur technique à la Société Télécommunications Radioélectriques et 
Téléphoniques pour avoir guidé ses travaux avec une grande clairvoyance et les 
avoir efficacement stimulés pendant de nombreuses années. L'auteur adresse aussi 
ses vifs remerciements à l'ensemble de ses collaborateurs pour leurs contributions 
et pour l’assistance constante qu'ils ont apportée. 
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INTRODUCTION 



Le signal est le support de l'information émise par une source et destinée à un 
récepteur; c’est le véhicule de l'intelligence dans les systèmes. Il transporte les 
ordres dans les équipements de contrôle et de télécommande, il achemine sur les 
réseaux l'information, la parole ou l'image. Il est particulièrement fragile et doit être 
manipulé avec beaucoup de soins. Le traitement qu'il subit a pour but d'extraire des 
informations, de modifier le message qu'il transporte ou de l'adapter aux moyens de 
transmission; c'est là qu'interviennent les techniques numériques. En effet, si l'on 
imagine de substituer au signal un ensemble de nombres qui représentent sa gran- 
deur ou amplitude à des instants convenablement choisis, le traitement, même dans 
sa forme la plus élaborée, se ramène à une séquence d’opérations logiques et arith- 
métiques sur cet ensemble de nombres, associées à des mises en mémoire. 

La conversion du signal continu analogique en un signal numérique est réali- 
sée par des capteurs qui opèrent sur des enregistrements ou directement dans les 
équipements qui produisent ou reçoivent le signal. Les opérations qui suivent cette 
conversion sont réalisées par des calculateurs numériques agencés ou programmés 
pour effectuer l’enchaînement des opérations définissant le traitement. 

Avant d'introduire le contenu des différents chapitres du présent ouvrage, il 
convient de donner une définition précise du traitement considéré. 

Le traitement numérique du signal désigne l'ensemble des opérations, calculs 
arithmétiques et manipulations de nombres, qui sont effectués sur un signal à traiter, 
représenté par une suite ou un ensemble de nombres, en vue de fournir une autre 
suite ou un autre ensemble de nombres, qui représentent le signal traité. Les fonctions 
les plus variées sont réalisables de cette manière, comme l'analyse spectrale, le filtrage 
linéaire ou non linéaire, le transcodage, la modulation, la détection, l'estimation et 
l'extraction de paramètres. Les machines utilisées sont des calculateurs numériques. 

Les systèmes correspondant à ce traitement obéissent aux lois des systèmes 
discrets. Les nombres sur lesquels il porte peuvent dans certains cas être issus d'un 
processus discret. Cependant, ils représentent souvent l'amplitude des échantillons 
d'un signal continu et dans ce cas, le calculateur prend place derrière un dispositif 
convertisseur analogique-numérique et éventuellement devant un convertisseur 
numérique-analogique. Dans la conception de tels systèmes et l’étude de leur fonc- 
tionnement, la numérisation du signal revêt une importance fondamentale et les 
opérations d'échantillonnage et de codage doivent être analysées dans leur prin- 
cipe et leurs conséquences. La théorie des distributions constitue une approche 
concise, simple et efficace pour cette analyse. Après un certain nombre de rappels 
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sur l'analyse de Fourier, les distributions et la représentation des signaux, le cha- 
pitre premier rassemble les résultats les plus importants et les plus utiles sur 
l'échantillonnage et le codage d'un signal. 

L'essor du traitement numérique date de la découverte d'algorithmes de calcul 
rapide de la Transformée de Fourier Discrète. En effet, cette transformation est à 
la base de l’étude des systèmes discrets et elle constitue dans ce domaine numé- 
rique l’équivalent de la Transformation de Fourier dans le domaine analogique, 
c'est le moyen de passage de l'espace des temps discret à l'espace des fréquences 
discret. Elle s'introduit naturellement dans une analyse spectrale avec un pas de 
fréquence diviseur de la fréquence d'échantillonnage des signaux à analyser. 

Les algorithmes de calcul rapide apportent des gains tels qu'ils permettent de 
faire les opérations en temps réel dans de nombreuses applications pourvu que 
certaines conditions élémentaires soient remplies. Ainsi, la Transformation de 
Fourier Discrète constitue non seulement un outil de base dans la détermination 
des caractéristiques du traitement et dans l'étude de ses incidences sur le signal, 
mais de plus, elle donne lieu à la réalisation d’équipements toutes les fois qu'une 
analyse de spectre intervient, par exemple, dans les systèmes comportant des bancs 
de filtres ou quand, par la puissance de ses algorithmes, elle conduit à une 
approche avantageuse pour un circuit de filtrage. Les chapitres 2 et 3 lui sont 
consacrés ; ils donnent d'une part une présentation des propriétés élémentaires et 
du mécanisme des algorithmes de calcul rapide et de leurs applications, et d'autre 
part, un ensemble de variantes associées aux situations pratiques. En tant que sys- 
tème, le calculateur de Transformée de Fourier Discrète est un système linéaire 
discret, invariant dans le temps. 

Une grande partie du présent ouvrage est consacrée à l'étude des systèmes 
linéaires discrets invariants dans le temps à une dimension, qui sont facilement 
accessibles et très utiles. Les systèmes à plusieurs dimensions et en particulier à 
deux et trois dimensions connaissent un grand développement; ils sont appliqués 
par exemple aux images; cependant, leurs propriétés se déduisent en général de 
celles des systèmes à une dimension dont ils ne sont souvent que des extensions 
simplifiées. Les systèmes non linéaires ou variables dans le temps, soit contiennent 
un sous-ensemble important qui présente les propriétés de linéarité et invariance 
temporelle, soit peuvent s'analyser avec les mêmes techniques que les systèmes 
ayant ces propriétés. 

La linéarité et l'invariance temporelle entraînent l'existence d'une relation de 
convolution qui régit le fonctionnement du système, ou filtre, ayant ces propriétés. 
Cette relation de convolution est définie à partir de la réponse du système au 
signal élémentaire que représente une impulsion, la réponse impulsionnelle, par 
une intégrale dans le cas des signaux analogiques. Ainsi, si x (t) désigne le signal à 
filtrer, h(t) la réponse impulsionnelle du filtre, le signal filtré y ( t) est donné par 
l’équation : 



y(t) = 



h (x) v (t - x) dx 
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Dans ces conditions, une telle relation qui pourtant traduit directement le fonc- 
tionnement réel du filtre, offre un intérêt pratique limité. En effet, d'une part il 
n’est pas très aisé de déterminer la réponse impulsionnelle à partir des critères qui 
définissent l'opération de filtrage envisagée et d'autre part une équation compor- 
tant une intégrale ne permet pas facilement de reconnaître et vérifier le comporte- 
ment du filtre. La conception est beaucoup plus facile à aborder dans le domaine 
des fréquences car la transformation de Laplace ou la transformation de Fourier 
permettent d'accéder à un plan transformé où les relations de convolution du plan 
amplitude-temps deviennent de simples produits de fonctions. A la réponse impul- 
sionnelle, la transformation de Fourier fait correspondre la réponse en fréquence 
du système, et le filtrage se ramène au produit de cette réponse en fréquence par la 
transformée de Fourier, ou spectre, du signal à filtrer. 

Dans les systèmes numériques, qui sont du type discret, la convolution se tra- 
duit par une sommation. Le filtre est défini par une suite de nombres qui constitue 
sa réponse impulsionnelle. Ainsi, si la suite à filtrer s’écrit x(n ), la suite filtrée y (n) 
s’exprime par la sommation suivante, où n et m sont des entiers : 

y ( n ) = S h (m)x ( n - m ) 

m 



Deux cas se présentent alors. Ou bien la sommation porte sur un nombre fini 
de termes, c’est-à-dire que les h ( m ) sont nuis sauf pour un nombre fini de valeurs 
de la variable entière m. Le filtre est dit à réponse impulsionnelle finie ; en faisant 
allusion à sa réalisation, on le désigne encore par non récursif car il ne nécessite 
pas de boucle de réaction de la sortie sur l'entrée dans sa mise en œuvre. Il est à 
mémoire finie, puisqu'il ne garde le souvenir d'un signal élémentaire, une impul- 
sion par exemple, que pendant une durée limitée. Les nombres h (m) sont appelés 
les coefficients du filtre, qu'ils définissent complètement. Ils peuvent se calculer 
d'une manière directe très simple, par exemple en faisant le développement en 
série de Fourier de la réponse en fréquence à réaliser. Ce type de filtre présente 
des caractéristiques originales très intéressantes; par exemple, la possibilité d'une 
réponse rigoureusement linéaire en phase, c’est-à-dire d'un temps de propagation 
de groupe constant; les signaux dont les composantes se trouvent dans la bande 
passante du filtre ne sont pas déformés à la traversée de ce filtre. Cette possibilité 
est exploitée dans les systèmes de transmission de données ou en analyse spectrale 
par exemple. 

Ou bien la sommation porte sur un nombre infini de termes, les h ( m ) ont une 
infinité de valeurs non nulles; le filtre est dit à réponse impulsionnelle infinie ou 
encore de type récursif, car il faut réaliser sa mémoire par une boucle de réaction 
de la sortie sur l'entrée. Son fonctionnement est régi par une équation selon 
laquelle un élément de la suite de sortie y (n) est calculée par la sommation pondé- 
rée d'un certain nombre d'éléments de la suite d'entrée x(n) et d'un certain 
nombre d'éléments de la suite de sortie précédents. Par exemple, si L et K sont des 
entiers, le fonctionnement du filtre peut être défini par l'équation suivante : 
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L K 

y(ri)= X a,[X (n - 1) - X b k y(n-k ) 

/ = 0 A = 1 



Les a ; (Z = 0, 1, L) et b k (k = 1, 2, .... K) sont les coefficients. Comme pour les 
filtres analogiques, l'étude de ce type de filtre ne se fait pas en général simplement 
de manière directe ; il est nécessaire de passer par un plan transformé. La transfor- 
mation de Laplace ou la transformation de Fourier pourraient être utilisées. 
Cependant, il existe une transformation beaucoup mieux adaptée, la transforma- 
tion en Z, qui est l'équivalent pour les systèmes discrets. Un filtre est caractérisé 
par sa fonction de transfert en Z, désignée généralement par H (Z), et qui fait 
intervenir les coefficients par l’équation suivante : 



H (Z) = 



L 



Z 

/ = 0 



a{L~ l 



K 

1 + X b k Z-* 

k = t 



Pour obtenir la réponse en fréquence du filtre, il suffit de remplacer dans H (Z) la 
variable Z par l'expression suivante où /désigne la variable fréquence et T la 
période d'échantillonnage des signaux : 

Z = e> 2% f' v 

Dans cette opération, à l'axe imaginaire, dans le plan de Laplace, correspond 
le cercle de rayon unité centré à l'origine dans le plan de la variable Z. Il apparaît 
clairement que la réponse en fréquence du filtre défini par H (Z) est une fonction 
périodique ayant pour période la fréquence d’échantillonnage. Une autre repré- 
sentation de la fonction H (Z) est utile pour la conception des filtres et l'étude d'un 
certain nombre de propriétés, celle qui fait apparaître les racines du numérateur 
appelées zéros du filtre, Z ; (/ = 1, 2, ..., L) et les racines du dénominateur appelées 
pôles, P*(fc = 1,2,..., K) : 

L 

n ri-Z/Z- 1 ) 

U{Z) = a 0 ^ 

n (î-p.z 1 ) 

Le terme a 0 est un facteur d’échelle qui définit le gain du filtre. La condition de sta- 
bilité du filtre s’exprime très simplement par la contrainte suivante : tous les pôles 
doivent être à l'intérieur du cercle unité. La position des pôles et des zéros par rap- 
port au cercle unité, permet une appréciation très simple et très utilisée des carac- 
téristiques du filtre. 

Un ensemble de quatre chapitres est consacré à l'étude des caractéristiques 
de ces filtres numériques. Le chapitre IV présente les propriétés des systèmes 
linéaires discrets invariants dans le temps, rappelle les propriétés principales de la 
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transformation en Z et donne les éléments nécessaires à l'étude des filtres. Le cha- 
pitre V traite des filtres à réponse impulsionnelle finie : leurs propriétés sont étu- 
diées, les techniques de calcul des coefficients sont décrites ainsi que les structures 
de réalisation. Les filtres à réponse impulsionnelle infinie étant généralement réali- 
sés par une mise en cascade de cellules élémentaires du premier et second ordre ; le 
chapitre VI décrit ces cellules et leurs propriétés, ce qui d'une part facilite considé- 
rablement l'approche de ce type de système et d'autre part fournit un ensemble de 
résultats très utiles dans la pratique. Le chapitre VII donne les méthodes de calcul 
des coefficients des filtres à réponse impulsionnelle infinie et traite les problèmes 
apportés par la réalisation, avec les limitations qu'elle implique et leurs consé- 
quences, en particulier le bruit de calcul. 

Les filtres à réponse impulsionnelle infinie ayant des propriétés comparables à 
celles des filtres analogiques continus, il est naturel d'envisager pour leur réalisa- 
tion des structures du même type que celles qui sont couramment employées en fil- 
trage analogique. C'est l'objet du chapitre VIII qui présente des structures en 
chaîne. Une digression est faite avec les dispositifs à commutation de capacités, qui 
ne sont pas de type numérique au sens strict, mais qui sont néanmoins de type 
échantillonné et sont des compléments très utiles aux filtres numériques. Pour gui- 
der l'utilisateur, un résumé des mérites respectifs des structures décrites est donné 
en fin de chapitre. 

Certains équipements, par exemple en instrumentation ou dans le domaine des 
télécommunications, font intervenir des signaux représentés par une suite de 
nombres complexes. Dans l'ensemble des signaux de ce type, une catégorie présente 
un intérêt pratique notable, celle des signaux analytiques. Leurs propriétés sont étu- 
diées au chapitre IX, ainsi que la conception des dispositifs adaptés à la génération 
ou au traitement de tels signaux. Des notions complémentaires sur le filtrage sont 
également données dans ce chapitre, qui présente, d'une manière unifiée, les princi- 
pales techniques d'interpolation. 

Les machines de traitement numérique, quand elles fonctionnent en temps 
réel, opèrent à une cadence étroitement liée à la fréquence d'échantillonnage des 
signaux et leur complexité dépend du volume d'opérations à faire et de l'intervalle 
de temps disponible pour les réaliser. La fréquence d’échantillonnage des signaux 
est généralement imposée à l’entrée ou à la sortie des systèmes, mais à l'intérieur 
du système lui-même, il est possible de la faire varier pour l'adapter aux caractéris- 
tiques du signal et du traitement, et ainsi de réduire le volume d'opérations et la 
cadence des calculs. Une simplification des machines, qui peut être très importante, 
est obtenue en adaptant tout au long du traitement la fréquence d'échantillonnage 
à la largeur de bande du signal utile, c’est le filtrage multicadence présenté au cha- 
pitre X. Les incidences sur les caractéristiques du traitement sont décrites ainsi que 
les méthodes de réalisation. Des règles d'utilisation et d'évaluation sont fournies. 
Cette technique donne des résultats particulièrement intéressants pour les filtres à 
bande passante étroite ou pour la mise en œuvre d'ensembles appelés bancs de 
filtres. Dans ce dernier cas, le système associe à un ensemble de circuits déphaseurs 
un calculateur de Transformée de Fourier Discrète. 
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Les bancs de filtres pour la décomposition et la reconstruction des signaux 
sont devenus un outil de base pour la compression. Leur fonctionnement est décrit 
aux chapitres 11 et 12 avec les méthodes de calcul et les structures de réalisation. 

Les filtres peuvent être déterminés à partir de spécifications dans le temps; 
c’est le cas par exemple de la modélisation d'un système comme décrit au chapitre 13. 
Si les caractéristiques varient, il peut être intéressant de modifier les coefficients en 
fonction des évolutions du système. Cette modification peut dépendre d'un critère 
d'approximation et se faire à une cadence qui peut atteindre la cadence d’échan- 
tillonnage du système ; alors le filtre est dit adaptatif. Le chapitre 14 est consacré 
au filtrage adaptatif, dans le cas le plus simple, mais aussi le plus courant et le plus 
utile, celui où le critère d'approximation retenu est la minimisation de l'erreur qua- 
dratique moyenne et où les variations des coefficients se font suivant l'algorithme 
du gradient. Après un ensemble de rappels donnés au chapitre 13 sur les signaux 
aléatoires et leurs propriétés, en particulier la fonction et la matrice d'autocorréla- 
tion dont les valeurs propres jouent un rôle important, l'algorithme du gradient est 
présenté au chapitre 14 et ses conditions de convergence sont étudiées. Ensuite les 
deux paramètres d'adaptation principaux, la constante de temps et l'erreur rési- 
duelle sont analysés, ainsi que la complexité arithmétique. Différentes structures 
de réalisation sont proposées. 

Pour terminer, le chapitre 15 décrit brièvement quelques applications, en 
montrant comment les méthodes et techniques de base sont exploitées. 
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Chapitre 1 



La numérisation du signal 
Échantillonnage et codage 



La conversion d'un signal analogique sous forme numérique implique une double 
approximation. D'une part, dans l’espace des temps, le signal fonction du temps 
s ( t ) est remplacé par ses valeurs 5 («T ) à des instants multiples entiers d'une durée 
T ; c’est l'opération d'échantillonnage. D'autre part, dans l'espace des amplitudes, 
chaque valeur 5 (zzT ) est approchée par un multiple entier d'une quantité élémen- 
taire q ; c’est l'opération de quantification. La valeur approchée ainsi obtenue est 
ensuite associée à un nombre ; c’est le codage, ce terme étant souvent utilisé pour 
désigner l’ensemble, c’est-à-dire le passage de la valeur ,s(nT) au nombre qui la 
représente. 

L'objet du présent chapitre est d'analyser l'incidence sur le signal de ces deux 
approximations. 

Pour mener à bien cette tâche, on utilise deux outils de base qui sont l'analyse 
de Fourier et la théorie des distributions. 



1.1 L'ANALYSE DE FOURIER 

L’analyse de Fourier est un moyen de décomposer un signal en une somme de 
signaux élémentaires particuliers, qui ont la propriété d’êtres faciles à mettre en 
œuvre et à observer. L’intérêt de cette décomposition réside dans la fait que la 
réponse au signal d’un système obéissant au principe de superposition peut être 
déduite de la réponse aux signaux élémentaires. Ces signaux élémentaires sont 
périodiques et complexes, afin de permettre une étude en amplitude et en phase 
des systèmes ; ils s’expriment par la fonction s e ( t ) telle que : 

s e (t) = ei lK f‘ = cos (2k ft) + j sin (2k ft) (1.1) 

où /représente l’inverse de la période, c’est la fréquence du signal élémentaire. 
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Dans la mesure où les signaux élémentaires sont périodiques, il est clair que 
l'analyse se simplifie dans le cas où le signal est lui-même périodique. Ce cas va 
être examiné d'abord, bien qu'il ne corresponde pas aux signaux les plus intéres- 
sants, puisqu'un signal périodique est parfaitement déterminé et ne porte pratique- 
ment pas d'information. 



1.1.1 Développement en série de Fourier d'une fonction périodique 

Soit s(t), une fonction de la variable Apériodique et de période T, c’est-à-dire satis- 
faisant la relation : 



s (A + T ) = s (A) 



( 1 . 2 ) 



Sous certaines conditions, on démontre que cette fonction est développable en 
série de Fourier, c'est-à-dire que l'égalité suivante est vérifiée : 



s (A) = X C n ei 2nntlT 



(1.3) 



L'indice n est un entier et les C n sont appelés les coefficients de Fourier; ils 
sont définis par l'expression : 



1 

r = — 

v -'n j 



s (t) e~> 2mtrr dt 



(1.4) 



En fait les coefficients de Fourier minimisent l'écart quadratique entre la fonc- 
tion s (t) et le développement (1.3). En effet la valeur (1.4) est obtenue en dérivant 
par rapport au coefficient d'indice n l'expression : 

(s(t)~ £ C m ei lmn,n \ dt 

. 0 \ m = -oo / 

et en annulant cette dérivée. 

Exemple : développement en série de Fourier de la fonction i p (t) constituée 
par une suite d'impulsions, séparées par la durée T, de largeur t et d'amplitude a, 
centrée sur l'origine des temps (fig. 1.1). 



ip(t) 




Fig. 1.1. Suite d'impulsions 
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Les coefficients C„ s’écrivent : 



C„ 



T 



sm nn — 
_ ai V T 

ae -]2nntlT fa = — - 

T T 



- t/2 



nn — 
T 



et le développement est donné par : 



sm \nn — 
CT * l T 

*1.(0= T £ - T 

i n = — oo L 

nn — 
T 



gj2nnt/T 



(1.5) 



(1.6) 



On imagine l'importance que prend cet exemple dans l’étude des systèmes 
échantillonnés. 

Les propriétés des développements en série de Fourier sont présentées dans l'ou- 
vrage [1]. Une propriété importante est exprimée par l’égalité de Bessel-Parseval 
qui traduit le fait que dans la décomposition du signal il y a conservation de la 
puissance : 



oo 



Z icj 2 



î 

T 



'T 
. 0 



s(/)| 2 dt 



(1.7) 



Les signaux élémentaires qui résultent de la décomposition d'un signal périodique 

1 

ont des fréquences qui sont des multiples entiers de —, l'inverse de la période; ils 

couvrent un ensemble discret de l'espace des fréquences. Par contre si le signal 
n’est pas périodique, les signaux élémentaires résultant de la décomposition cou- 
vrent un domaine continu de l'espace des fréquences. 



1.1.2 Transformation de Fourier d'une fonction 

Soit s(t) une fonction de la variable f ; sous certaines conditions on démontre l’éga- 
lité suivante : 



avec 



s(t) = 



S (f)eWdf 



(1.8) 



S (/) = 



s(t) e~l 2K fi dt 



(1.9) 



La fonction S (/) est la transformée de Fourier de s(t). Plus communément S (/) est 
appelé spectre du signal s (f). 
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Exemple : soit à calculer la transformée de Fourier I(/) d'une impulsion isolée 
i(t) de largeur T, d'amplitude a et centrée sur l'origine des temps (fig. 1.2) 



!(/) = 



i ( t ) e~i 2n f‘ dt= a 



e-J 2n ft dt 



J-l/2 



I(/) = UT 



sin (ît/î) 
nfz 



( 1 . 10 ) 



Fig. 1.2. Impulsion isolée 




La figure 1.3 représente la fonction I(/), qui sera très fréquemment utilisée 
par la suite. Il est important de remarquer qu’elle s'annule aux fréquences mul- 
tiples entiers non nuis de l'inverse de la durée de l'impulsion. 

L'Annexe 1 donne une tabulation de cette fonction. 

La correspondance entre coefficients de Fourier et spectre apparaît nettement 
sur cet exemple. En effet, en rapprochant les relations (1.6) et (1.10) on vérifie que, 
1 

au facteur — près, les coefficients du développement en série de Fourier de la 

suite d'impulsions correspondent aux valeurs que prend le spectre de l'impulsion 
isolée aux fréquences multiples entiers de l'inverse de la période des impulsions. 

En fait, on a la relation : 

1 ! n 

C " = T S ( T 

Une relation comparable à l'égalité de Bessel-Parseval existe pour une fonc- 
tion non périodique. Dans ce cas, c’est non plus la puissance mais l’énergie du 
signal qui se trouve conservée : 



|s(/)M/= 



s (f)| 2 dt 



( 1 . 11 ) 



Soit s\t) la dérivée de la fonction 5 ( t ) ; sa transformée de Fourier S d {f) s’écrit : 



S Af) 



e~i 2n ft .s\t) dt- j2nf.S(f) 



( 1 . 12 ) 



Ainsi prendre la dérivée d'un signal amène une multiplication de son spectre par 
7'2tc/. 

Une propriété essentielle de la transformation de Fourier, qui est en fait la 
principale raison de son utilisation, est qu’elle transforme un produit de convolu- 
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tion en un produit simple. En effet soit deux fonctions du temps x(t ) et h (t) dont 
les transformées de Fourier sont respectivement X(/) et H(/). Le produit de 
convolution y (/) est défini par : 



y(t) = x{t)*h{t) = 



x(t-x)h(x) dx 



(1.13) 



La transformée de Fourier de ce produit s’écrit : 



Y (/) = 



x (t - t) h (t) di e~i 2K ft dt 



Y (/) = 



h(x) e~i 2n f x dx. 



x(u) e~l 2n f u du = H (/) .X(/) 



Réciproquement, on montre que la transformée de Fourier d'un produit simple est 
un produit de convolution. 

Un résultat intéressant pour l'étude de l’échantillonnage et se rapportant à 
l'exemple ci-dessus peut être déduit directement de ces propriétés. En effet soit à 
calculer la transformée de Fourier II(/) de la fonction i 2 (t) ; d'après les relations 
(1.10) et (1.13), il vient : 

II(/)=I(/)*I(/) = «.I(/) (1.14) 

et par suite : 

sin (rapt) sin [jc(/-<p) t] 1 sin (nfx) 

. tupi ‘ 7t(/-(p)T ^ ~ x n fx 

n 

En prenant /= - , pour tout entier n non nul, on a : 

f" sin (rapt) sin [jt (cpr - n)] 

. r — dm - 0 (1.15) 

7tcpx 7t(cpx-n) Y v ' 



sin n (x - n) 

Les fonctions ; — , avec n entier, forment un ensemble de fonctions ortho- 

x! (x - n) 

gonales. 

La définition et les propriétés de la transformation de Fourier s’étendent aux 
fonctions de plusieurs variables. Soit s(x x , x 2 , ..., x n ) une fonction de n variables 
réelles, la transformée de Fourier est une fonction S(X l5 X 2 , ..., définie par : 

s(^ 2 ,..a„) 




s(x v x 2 , ...,x n )e~ j2n (\ x i + x 2*2 + — + K x n> dx 1 dx 2 , ..., dx n (1.16) 
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Si la fonction s (xq, x 2 , ..., x n ) est séparable c’est-à-dire si : 
s{x l ,x 2 ,..., x„ ) = s (x 2 ) s (x 2 ) ... s (x n ) alors il vient : 

S(^A 2 , —,K) = S(^i) S(^)...S (À„) 

Les variables x, (1 i =£ n) représentent souvent des distances, par exemple 
dans le cas bidimensionnel, et les X, sont alors appelées fréquences spatiales. 

Dans l'étude des signaux échantillonnés, la transformation de Fourier va être 
appliquée aux distributions. 



1.2 LES DISTRIBUTIONS 

Les distributions mathématiques constituent une définition mathématique correcte 
des distributions rencontrées en physique [1]. 

1.2.1 Définition 

On appelle distribution D une fonctionnelle linéaire continue sur l'espace vectoriel 
2) des fonctions définies sur IR", indéfiniment dérivables et à support borné. 

A toute fonction cp appartenant à 2), la distribution D associe un nombre com- 
plexe D((p), qui sera aussi noté par (D, cp>, avec les propriétés : 

- D( 9l + (p 2 ) = + D(cp 2 ). 

- D(^cp) = ?iD(cp) où X est un scalaire. 

- Si <p ( converge vers cp quand j tend vers l'infini, la suite D ( cp ■) converge vers 

D(q>). 

Exemples : 

• Si f(t ) est une fonction sommable sur tout ensemble borné, elle définit une 
distribution par : 

(D / ,<P)= /(O <p(f) dt (1.17) 

• Si cp' désigne la dérivée de (p, la fonctionnelle : 

(D, <p) = dt = (f,q') 

est une distribution. 

• La distribution de Dirac 5 est définie par : 

<ô, cp) = <p (0) 

La distribution de Dirac au point réel x est définie par : 

<ô(f-x),cp) = cp(x) (1.20) 

On dit que cette distribution représente la masse + 1 au point x. 



(1.18) 

(1.19) 




© Dunod. La photocopie non autorisée est un délit. 



1 .2 Les distributions 



13 



• Soit l'impulsion i(t) de durée T, d'amplitude a = 1/x, centrée sur l'origine. 
Elle définit une distribution D, : 

1 t/ 2 

(D,-,<P>=- 9(0 àt 

1 J-x/2 

Pour des valeurs de x très petites on obtient : 

(D;, 9) — 9(0) 

c'est-à-dire que la distribution de Dirac peut être considérée comme la limite, 
quand x tend vers 0, de la distribution D, . 

1.2.2 Dérivation des distributions 

dD 

On définit la dérivée — d'une distribution D par la relation : 
dt 

, <3D , 09 

<^>=-< D -i> < 121 » 

Soit par exemple la fonction Y de Heaviside, ou échelon unité, égale à 0 si 
t < 0 et + 1 si t 5* 0. 

<^,9> = -<Y, -^> = -[ 9' (f) df = 9 (0) = <ô, 9) (1.22) 

ot ot J 0 

Il en résulte que la discontinuité de Y apparaît sous la forme d'une masse ponc- 
tuelle unitaire dans sa dérivée. 

Cet exemple illustre un intérêt pratique considérable de la notion de distribution, 
qui permet d'étendre aux fonctions discontinues un certain nombre de concepts et 
de propriétés des fonctions continues. 

1.2.3 Transformation de Fourier d'une distribution 

Par définition la transformée de Fourier d'une distribution D est une distribution 
notée FD telle que : 

(FD, 9) = (D, F9) (1.23) 

Par application de cette définition aux distributions à support ponctuel il vient : 

<FÔ, 9) = <5, F9> = [ 9 (0 dt = { 1, 9) (1.24) 

Par suite : F8 = 1. 

De même F8 ( t-a ) = e~i 2n f a . 
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Un cas fondamental pour l'étude de l'échantillonnage est celui que constitue 
la suite des distributions de Dirac décalées de T, notée u et telle que : 

oo 

u (t) = X ô(f-nT) (1.25) 

n = - oo 

Cette suite est une distribution de masses unitaires aux points dont l'abscisse est un 
multiple entier de T. Sa transformée de Fourier s’écrit : 

oo 

F u= X e~> 2% f nT = U(/) (1.26) 

n = - oo 

On démontre que cette somme est en fait une distribution ponctuelle. 

Une démonstration intuitive peut être obtenue à partir du développement en 
série de Fourier de la fonction i p (t) constituée par la suite d'impulsions séparées 
par la durée T, de largeur x, d'amplitude 1/x, centrée sur l'origine des temps. 

En effet on peut considérer que : u (t) = lim i (t). 



En se reportant à la relation (1.6) on trouve : 




oo 

g-j2nnt/T 

n = -oo 



Il en résulte que : 



00 1 °° ( n \ 

u (/) = X e~’W=- 2 s/-- 

n = - oo 1 n = — co \ 1/ 



(1.27) 



Cette propriété fondamentale démontrée dans l'ouvrage [1], ainsi que dans l'ou- 
vrage [2], s’exprime comme suit : 

La transformée de Fourier de la distribution temporelle comportant une 
masse unitaire en chaque point dont l’abscisse est un multiple entier de T est une 
distribution fréquentielle comportant la masse 1/T aux points dont l’abscisse est un 
multiple entier de 1/T. 

Ce résultat va être utilisé pour étudier l’échantillonnage d'un signal. 

La propriété que possède la transformation de Fourier d'échanger convolu- 
tion et multiplication s’applique également aux distributions. 

Avant d'étudier les incidences sur le signal des opérations d'échantillonnage 
et quantification, il est utile de caractériser les signaux qui sont les plus fréquem- 
ment traités. 



1 .3 LES PRINCIPAUX SIGNAUX TRAITÉS 



Les signaux sont définis par une fonction du temps s(t). Cette fonction peut être 
une expression analytique ou la solution d'une équation différentielle, auquel cas 
le signal est appelé déterministe. 
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1 .3.1 Les signaux déterministes 

Les signaux de ce type les plus utilisés sont les signaux sinusoïdaux ; par exemple : 

s(t) - A cos (cot + a) 

où A est l'amplitude, co = 2ît/la pulsation et a la phase du signal. 

Ils sont faciles à reproduire, à reconnaître aux différents points d'un système 
et offrent une possibilité de visualisation simple des caractéristiques. De plus, 
comme indiqué aux paragraphes précédents, ils servent de base à la décomposition 
d'un signal déterministe quelconque, par l'intermédiaire de la Transformation de 
Fourier. 

Si le système considéré est linéaire et invariant dans le temps, il peut être carac- 
térisé par sa réponse en fréquence FI(cû). Pour chaque valeur de la fréquence, H(co) 
est un nombre complexe dont le module est l'amplitude de la réponse. Par conven- 
tion on désigne par phase de la réponse du système la fonction cp (co) telle que : 

H(co) = |H(co)| <rM®) (1.28) 

Cette convention permet d'exprimer le temps de propagation de groupe t(cû), 
fonction positive dans les systèmes réels, par : 

dm 

x(co) = -g (1.29) 

Le temps de propagation de groupe fait référence aux lignes de transmission, 
sur lesquelles les différentes fréquences d'un signal se propagent à des vitesses dif- 
férentes, ce qui entraîne une dispersion dans le temps de l'énergie du signal. Pour 
illustrer cette notion, soit deux fréquences proches co ± Aco auxquelles correspon- 
dent les phases par unité de longueur cp ± Acp. Le signal somme s’écrit : 

s (t) = cos [(co + Aco)t - (cp + Acp)] + cos [(co - Aco)t - (cp - Acp)] 



ou encore 

s(t) =2 cos (cof - cp) cos (Acot - Acp) 

C’est un signal modulé et il n'y a pas de dispersion si les deux facteurs subis- 
sent le même retard par unité de longueur, c’est-à-dire si Acp/ Aco est une constante. 
Le temps de propagation de groupe caractérise donc la dispersion apportée à un 
signal par une ligne de transmission ou un système équivalent. 

En appliquant au système le signal sinusoïdal s(t), on obtient en sortie le 
signal résultant s r ( t) tel que : 

s r (t) = A. |H(co)| cos [coï + a-cp(co)] (1.30) 

C’est encore un signal sinusoïdal et la comparaison avec le signal appliqué 
permet une visualisation de la réponse du système. On imagine aisément l'impor- 
tance de cette procédure pour les opérations de test par exemple. 

Les signaux déterministes cependant ne représentent pas très bien les signaux 
réels, car, en fait, ils ne portent pas d'information, si ce n’est pas leur présence 
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même. Les signaux réels sont généralement caractérisés par une fonction s(t ) aléa- 
toire. Pour le test et l'analyse des systèmes on utilise aussi des signaux aléatoires, 
mais qui présentent des caractéristiques particulières pour ne pas compliquer exa- 
gérément la génération et l'exploitation. Une étude des signaux aléatoires est faite 
dans le tome 2 de la référence [2]. 



1 .3.2 Les signaux aléatoires 



Un signal aléatoire est défini à chaque instant t par la loi de probabilité de son 
amplitude s(t). Cette loi peut s’exprimer par une densité de probabilité p(x, t ) 
définie comme suit : 



p (X t) = lim 

^ v ' A.V — > o 



Proba [x =£ s (t ) « x + Ax] 
Ax 



(1.31) 



Il est stationnaire si ces propriétés statistiques sont indépendantes du temps, 
c'est-à-dire que sa densité de probabilité est indépendante du temps : 

p{x, t) =p(x) 



Il est du second ordre s'il possède un moment d'ordre 1 appelé valeur 
moyenne, qui est l'espérance mathématique de s(t), notée E[s(f)] et définie 
par : 



m\ (t) = E[.y(f)] 



x.p ( x , t ) dx 



et un moment d'ordre 2, appelé fonction covariance : 



E[s(t ] ).s(t 2 )] = m 2 (t ] ,t 2 ) = 



x 1 .x 2 ■ p (x,, x 2 ; t], t 2 ) dx 1 dx 2 



où p(x 2 , x 2 ; q, t 2 ) est la densité de probabilité du couple de variables aléatoires 

[s(q),s(q)]. 

Le caractère de stationnarité peut être limité aux moments du premier et du 
second ordre ; on dit alors que le signal est stationnaire d'ordre 2 ou stationnaire au 
sens large, et pour un tel signal il vient : 



E(î(t)] 



x.p (x) dx = m l 



L'indépendance du temps se traduit comme suit pour le densité de probabilité 
p(X|,x 2 ; f : , t 2 ) : 

P (x b x 2 ; q, t 2 ) = P (x b x 2 ; 0, t 2 - q) =p(x,,x 2 ; x) 

avec 

x = q-q 

Seul intervient l'écart entre les deux instants d'observation du signal : 

E[(.s'(q)..s(q)] = m 2 (x) 



(1.33) 
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La fonction r xx (t) telle que : 

'■„W = E[î(0-s(Lt)] (1-34) 

prend le nom de fonction d'autocorrélation du signal aléatoire, qu’elle caractérise. 

Un signal aléatoire s(t) possède aussi une moyenne temporelle m T , qui est 
une variable aléatoire définie par : 

1 pm 

m T = lim — s(t) dt (1.35) 

T ^°° T -T/2 

L'ergodicité de cette moyenne exprime le fait qu'elle prend une valeur déterminée 
k avec la probabilité 1. Pour un signal stationnaire, l'ergodicité de la moyenne tem- 
porelle entraîne l’égalité avec la moyenne des amplitudes à un instant donné. En 
effet prenons l'espérance de la variable m T : 

r -[ rT/2 n ^ t/2 

E[m T ] = k=E lim — s(t) dt = lim — E [s(r)] dt = m 1 

T — > co ’X' _ t/2 T — ^ 00 ’T — T/2 

Ce résultat a des conséquences pratiques importantes puisqu'il fournit un moyen 
d’accéder aux propriétés statistiques du signal à un instant donné à partir de l’ob- 
servation de ce signal au cours du temps. 

L'ergodicité de la covariance dans le cas stationnaire est également très intéres- 
sante car elle conduit à la relation : 

i rT/2 

r xx (x) = lim s{t)s{t-x)dt (1.36) 

T_>0 ° T J-T/2 

La fonction d'autocorrélation du signal s(t ), r xx (t) est fondamentale pour l'étude 
des signaux stationnaires d'ordre deux ergodiques. Ses principales propriétés sont 
les suivantes : 

- C’est une fonction paire : 

- Son maximum est à l'origine et correspond à la puissance du signal P : 

r xx (0) = E[.C(0] = P 

- La densité spectrale de puissance est la tranformée de Fourier de la fonc- 
tion d'autocorrélation : 

&XX (/) = r xx (?) e~ i2nfl dt = 2 r xx (t) cos (2ji/c) dt 

. -■» o 

En effet : r xx (t) = j(t) * s(- t) et, si S(/) désigne la transformée de Fourier de s(t), 
il vient : 

*xx(f) = S(/).S(7) = |S(/)| 2 (1.37) 

Cette dernière propriété se traduit physiquement par le fait que plus le signal 
est à variation rapide, c’est-à-dire plus son spectre s’étend vers les fréquences éle- 
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vées, plus sa fonction d'autocorrélation est étroite. A la limite le signal est pure- 
ment aléatoire et la fonction s'annule pour T + 0. On se trouve en présence d'un 
signal appelé bruit blanc, et tel que : 

rjt) = Pô 

Alors la densité spectrale est constante : 

<M/) = p 

En fait un tel signal n'a pas de réalité physique puisque sa puissance est infinie, 
mais il constitue un modèle mathématique commode pour les signaux dont la den- 
sité spectrale est quasi constante sur une large bande de fréquence. 



1.3.3 Les signaux gaussiens 



Parmi les lois de probabilité que l'on peut considérer pour un signal .v ( t) . il est une 
catégorie qui présente un grand intérêt, celle des lois normales ou lois de Gauss. 
En effet les distributions aléatoires normales conservent leur caractère normal 
dans toute opération linéaire, par exemple la convolution par une distribution cer- 
taine, le filtrage, la dérivation ou l'intégration. Aussi ces distributions aléatoires 
sont-elles très utilisées pour la modélisation des signaux réels et le test des 
systèmes. 

Une variable aléatoire x est dite gaussienne si sa loi de probabilité a une den- 
sité p (x) qui suit la loi normale ou loi de Gauss : 

\ ( x - m ) 2 

pW= ?v5U~ (L38) 



La valeur m est la moyenne de la variable x; la variance o 2 est le moment d'ordre 
deux de la variable centrée (x - m) ; o est aussi appelé l'écart-type. 

( x — m\ 

La variable I I est dite réduite, elle a une moyenne nulle et un écart-type 

unité. Une tabulation et une représentation très utile sous forme de courbe sont 
fournies en annexe II. 

Une variable aléatoire est caractérisée par la loi de probabilité de son ampli- 
tude, mais aussi par l’ensemble de ses moments m n , tels que : 



m 



n 



x n p (x) dx 



(1.39) 



Ces moments sont les coefficients du développement en série entière d'une fonc- 
tion F(») appelée fonction caractéristique de la variable aléatoire x et définie par : 



F («) 




e> ux p (x) dx 



(1.40) 
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C’est, à un changement de variable près, la transformée de Fourier inverse de la 
densité de probabilité p (x) et l'on a également : 



p{x) = 



1 

2k 



e i“ x F(u) du 



(1.41) 



À partir de la relation (1.40) on obtient le développement en série entière sui- 
vant : 

0° (il An 

Hu) = Z { ~ m n (1.42) 

« = o n\ 

Et pour une variable gaussienne centrée : 

î 

F(m) = e~2° 2 " 2 (1.43) 



Par développement en série et identification avec (1.42), on obtient tous les 
moments : 



ru 



2 fl 



(2»)! 

77 ! 2" 



Par exemple, pour n - 2, on obtient m 4 - 3cri. Tous les moments d'ordre 
impair d'une variable gaussienne centrée sont nuis, d'après la définition de la loi de 
probabilité elle-même. 

La loi normale se généralise aux variables aléatoires à plusieurs dimensions 
[3]. La fonction caractéristique d'une variable gaussienne à k dimensions v (x, , x k ) 
s’écrit : 

1 k k 

F(«i, u k ) = e “2 jh r n u ‘ u i (1.44) 



avec : 

C = E(x,x ; ) 

La densité de probabilité est obtenue par transformation de Fourier. Dans le 
cas à 2 dimensions il vient : 



P (*b x 2 ) 



1 



2kg, G-, Vl - ‘ 



e 2 ( 1 -P) 



x\ 2rx l x 2 x\ 
o? o.o, 



(1.45) 



où r désigne le coefficient de corrélation : 

E (x, x 2 ) 



G, G 



1^2 



Un signal aléatoire .v ( t) est dit gaussien, si pour un ensemble de k instants 
/,( 1 ^ ^ k) la variable aléatoire à k dimensions s = [s(q), s(f fe )] est gaussienne. 

D'après la relation (1.44), la loi de probabilité de cette variable est complète- 
ment définie par la fonction d'auto-corrélation r xx (x) du signal s(t). 
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Exemple : 

Le signal défini par les équations suivantes : 

r xx (t) = a 2 e~RC 
1 * 



P (X) : 



a V2ît 



e 2o 2 



(1.46) 

(1.47) 



est une approximation d'un bruit gaussien blanc d'utilisation courante dans l'ana- 
lyse des systèmes ou la modélisation des signaux. C’est un signal stationnaire de 
moyenne nulle dont la densité spectrale n’est pas rigoureusement constante, mais 
correspond à une répartition uniforme filtrée par un filtre passe-bas de type RC. Il 
s’obtient par amplification du bruit d'agitation thermique aux bornes d'une résis- 
tance. 

La distribution normale peut être obtenue à partir d'une distribution de pro- 
babilité p(x) uniforme sur l'intervalle [0, 1], En effet soit p (y) la distribution dite 
de Rayleigh : 

y _zl 

P(y)= e 2 a 2 ; y 5*0 (1.48) 



qui a pour moment d'ordre deux ou puissance, 2o 2 , pour moyenne 
variance ^2 - ^ jet 2 . Par un changement de variable tel que : 

p (x) dx = p (y) dy 




a et pour 



il vient : 



d'où : 



p(y) = p (*) 



x = e 



dx dx 
dy dy ’ 

Zl 

2a 2 ; y = G 



dx 

dy 



y - — 

e 2 a 2 

a- 




(1.49) 



La distribution normale est obtenue en considérant deux variables y et x indé- 
pendantes et en posant : 

z — y cos 2îtx (1.50) 

La démonstration fait intervenir la variable : 

z' = y sin 2ttx 

En effet, en utilisant la correspondance entre coordonnées polaires et carté- 
siennes, on peut écrire : 

p (z, z') dz dz' =p(z) p (z') dz dz' = p(y) p (x) dx dy = p ( z ) p ( z ') ydy 2n dx 
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1.3.4 Facteur de crête d'un signal aléatoire 

Un signal aléatoire est défini à chaque instant par une loi de probabilité de son 
amplitude, souvent telle que cette amplitude n’est pas bornée. C’est le cas des 
signaux gaussiens, comme le montre la relation (1.38). 

Or le traitement d’un signal ne peut se réaliser que pour une gamme d’ampli- 
tudes limitée et des opérations de cadrage interviennent. Un paramètre important 
est le facteur de crête défini pour le signal comme le rapport d’une certaine ampli- 
tude A m à la valeur efficace a. Par convention cette amplitude A m est souvent prise 
comme la valeur qui n’est pas dépassée pendant plus de 10 5 du temps. Ce rapport 
est exprimé en décibels (dB) par F c tel que : 

F c = 201og[^] (1.51) 

où log désigne le logarithme en base 10. 

Pour un signal gaussien le facteur de crête est de 12,9 dB. Appliquée à un 
signal sinusoïdal cette définition conduit à un facteur de crête de 3 dB. 

Un modèle stationnaire utilisé pour représenter le signal téléphonique est 
constitué par le signal aléatoire dont la densité de probabilité des amplitudes suit 
la loi exponentielle, ou de Laplace, suivante : 

1 r M 

‘ (L52) 

Le facteur de crête dans ce cas s’élève à 17,8 dB. 

En conclusion, les fonctions aléatoires stationnaires d’ordre deux ergodiques, 
caractérisées par une loi de probabilité des amplitudes et une fonction d’autocorré- 
lation, permettent de modéliser la plupart des signaux à traiter et sont très utilisées 
dans l’étude et l’analyse des systèmes. 

En plus des possibilités de représentation des signaux il est important de pou- 
voir disposer d’une mesure globale, par exemple afin de pouvoir suivre un signal 
au cours du traitement. Une telle mesure est obtenue en définissant des normes sur 
la fonction qui représente le signal. 
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1.4 NORMES D'UNE FONCTION 

Une norme est une fonction positive réelle, qui vérifie les relations : 

|M| s* 0; fc||x|| = ||fcr|| 



où k est un réel positif. 

Une catégorie très utilisée de normes est l'ensemble des normes dites 
normes-Lp [4] : 

La norme-Lp d'une fonction continue s ( t) définie sur l'intervalle [0, 1] est 
notée ||s|| p et définie par : 



LJ 0 



i(t)\ p dt 



(1.53) 



Trois valeurs de p sont intéressantes : 

-P = l: 

i-i 

|.s(/)|r/r (1.53-fl) 

o 

- p = 2 : 

\s{t)\ 2 dt (1.53-6) 





c'est l'expression de l'énergie du signal s{t) 

- p = oo : 

IML= max KOI 



(1.53-c) 



Cette norme est aussi appelée norme de Tchebycheff. Les normes sont utilisées 
également dans les techniques d'approximation pour mesurer l'écart entre une 
fonction f (x) et la fonction à approcher F(x). L'approximation est faite au sens des 
moindres carrés si la norme L 2 est utilisée et au sens de Tchebycheff si la norme L. r 
est utilisée. 

Les normes -Lp peuvent être généralisées par l'introduction d'une fonction de 
pondération réelle positive p (x). La norme-Lp pondérée de la fonction d’écart 
f(x) - F(x) s’écrit alors : 



ll/(*)-F(*)L = 



\f(x)-F(x)\Pp(x) dx 



(1.53 -d) 



Ces notions sont appliquées dans le calcul des coefficients des filtres et aussi des 
facteurs d'échelle qui commandent les cadrages des données dans les mémoires. 
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1.5 L'OPÉRATION D'ÉCHANTILLONNAGE 

L'échantillonnage consiste à représenter un signal fonction du temps s(t) par ses 
valeurs ,v (zzT ) à des instants multiples entiers d'une durée T, appelée période 
d’échantillonnage. Une telle opération s’analyse de façon simple et concise par 
l'intermédiaire de la théorie des distributions. En effet, par définition, la distribu- 
tion de masses unitaires aux points de l'axe réel multiples entiers de la période T, 
associe à la fonction ,v ( t) l'ensemble de ses valeurs .v («T ) où n est un entier. 
Conformément aux notations précédemment retenues cette distribution est notée 
u ( t ) et s’écrit : 

oo 

u ( t ) = X ô(f — zzT ) 

n = - oo 

L'opération d'échantillonnage affecte le spectre S(/) du signal. Considérant la 
relation fondamentale (1.27), il apparaît que le spectre U(/) de la distribution u ( t) 

est constitué de raies d'amplitude — aux fréquences qui sont des multiples entiers 

— . Par suite u(t) s’exprime comme une 



X e> 2mtlT (1.54) 

n = - oo 

Alors la suite des valeurs de signal s (nT ) correspond au produit de l'ensemble 
des signaux élémentaires qui constituent u(t) par le signal s (t). C’est-à-dire que 
physiquement, l'opération d'échantillonnage est une modulation en amplitude par 
le signal d'une infinité de porteurs à des fréquences qui sont des multiples entiers 
de la fréquence d'échantillonnage f e = 1/T. Par suite le spectre du signal échan- 
tillonné comprend la fonction S (/), désignée par la bande de base, ainsi que les 
bandes images qui correspondent à la translation de la bande de base de multiples 
entiers de la fréquence d'échantillonnage. 

L'opération d'échantillonnage et son incidence sur le spectre du signal sont 
représentées sur la figure 1.4. 



de la fréquence d'échantillonnage fe = 

somme de signaux élémentaires : 

, . 1 
“(0 = Tf 



Fig. 1 . 4 . Incidence spectrale de l’échantillonnage 
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Le spectre du signal échantillonné S e (f) a pour expression le produit de 
convolution de S (/) par U(/) soit : 

w-ï-ifî) (l55) 



Il est important de remarquer que la fonction S e {f) est périodique, c’est-à-dire 
que l'échantillonnage a introduit une périodicité dans l’espace des fréquences, ce 
qui constitue une caractéristique fondamentale des signaux échantillonnés. 

L'opération d’échantillonnage telle qu’elle vient d’être décrite et que l'on 
désigne par échantillonnage idéal, peut sembler peu réaliste, dans la mesure où il 
apparaît difficile dans la réalité d'atteindre, de manipuler ou de restituer une valeur 
d'un signal à un instant ponctuel ; les échantillonnées réels ou les circuits qui resti- 
tuent les échantillons possèdent un certain temps d'ouverture. En fait on peut mon- 
trer que l'échantillonnage ou la restitution d'échantillons par des impulsions ayant 
une largeur donnée, introduit simplement une modification du spectre du signal. 

En effet dans l'opération d'échantillonnage du signal s(t) par la suite d'impul- 
sions séparées par la durée T, de largeur x et d'amplitude a, il se peut que l'on 
recueille à la période n une quantité a n qui s’écrit : 



’nT + t/2 
nT - x/2 



s (t) dt 



Cette quantité exprime le résultat de la convolution du signal s(t ) par l'impul- 
sion élémentaire i (t) et la fonction dont on prélève dans ce cas les valeurs aux ins- 
tants d'échantillonnage n T est la fonction s * i\ c’est-à-dire que le signal échan- 
tillonné a pour spectre non pas S (/) mais le produit : 



S (f)-ax. 



sin (jc/x) 
nfx 



Le raisonnement est le même pour le cas de la restitution d’échantillons avec 
une durée x. En fait c’est le produit de convolution des échantillons s(«T) avec 
l'impulsion élémentaire i (t) qui est restitué. 

D'où la proposition : 

L’échantillonnage ou la restitution d’échantillons par des impulsions de largeur 
t peut être traité comme un échantillonnage idéal ou une restitution idéale, à la 
condition de multiplier le spectre du signal par le spectre de l’impulsion élémentaire. 

En pratique dès que x est faible devant la période T la correction devient 
négligeable. 



1.6 L'ÉCHANTILLONNAGE EN FRÉQUENCE 

L’échantillonnage considéré ci-dessus est de type temporel. Cependant les proprié- 
tés énoncées sont aussi applicables à un échantillonnage de type fréquentiel. 
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Calculons le spectre d'une fonction périodique s p (t ) de période T. Une telle 
fonction peut être considérée comme résultant du produit de convolution de la 
fonction s(t), qui prend les valeurs de s p (t) sur une période et s'annule en dehors, 
et de la distribution ponctuelle u ( t) . Il en résulte la relation suivante entre les 
transformées de Fourier : 

S,(/) = U(/).S(/)= 1 n £ œ s(^)ô(/- ( 1 . 56 ) 

En fait on retrouve les coefficients du développement en série de Fourier de la 
fonction s p (t). Le cas où s(t) est une impulsion est représenté sur la figure 1.5. 

Il apparaît que le spectre de la fonction périodique s p ( t ) est un spectre de raies 
qui constituent un échantillonnage du spectre de la fonction prise sur une période. 
L'échantillonnage dans l'espace des fréquences correspond à une périodicité dans 
l'espace des temps. Cette interprétation est utile dans l'analyse numérique des 
spectres. 




T 







t 




1.7 LE THÉORÈME DE L'ÉCHANTILLONNAGE 

Ce théorème exprime les conditions dans lesquelles la suite des échantillons d'un 
signal représente correctement ce signal. Un signal est supposé être correctement 
représenté par la suite de ses échantillons prélevés avec la périodicité T s'il est pos- 
sible, à partir de cette suite de valeurs, de restituer intégralement le signal 
d'origine. 

L'échantillonnage a introduit une périodicité du spectre dans l'espace des fré- 
quences; restituer le signal d'origine, c’est supprimer cette périodicité, c’est-à-dire 
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éliminer les bandes images, opération qui peut être réalisée à l'aide d'un filtre 
passe bas dont la fonction de transfert H(/) vaut l/f e jusqu’à la fréquence fJ2 et 0 
aux fréquences supérieures. En sortie d'un tel filtre apparaît un signal continu, 
qu'il est possible d'exprimer en fonction des valeurs ^ (/zT ). La réponse impulsion- 
nelle du filtre h ( t ) s’écrit, d'après la relation (1.10) : 



MO = 



sin {ntl T) 
nt/T 



Le signal de sortie du filtre, s(t), correspond au produit de convolution de la 
suite 5 («T ) par la fonction h (t), soit : 



MO = 



X s(0)S(0-uT) 

n = - oo 



sin n(t- 0)/T 
Jc(f-0)/T 



dQ 



d'où : 



00 

s (0 = X s (nT ) 

n = - oo 



sin jc(f/T - n ) 
n(t/T - n ) 



(1.57) 



C’est la formule de calcul des valeurs du signal aux instants situés entre les 
échantillons. Pour les multiples de la période T elle fournit bien s («T). Le proces- 
sus de reconstitution du signal est représentée sur la figure 1.6. 




Fig. 1.6. Reconstitution du signal après échantillonnage 

Pour que le signal calculé s(t) soit identique au signal d'origine, il faut que le 
spectre S (/) soit identique au spectre du signal d'origine. Comme le montre la 
figure 1.6 cette condition est vérifiée si et seulement si le spectre d'origine ne 
contient pas de composantes aux fréquences supérieures ou égales à f e /2. 

Si ce n’est pas le cas, les bandes images chevauchent la bande de base comme 
sur la figure 1.7, on dit qu'il y a repliement de bande, et le filtre de restitution four- 
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nit un signal différent du signal d'origine. D'où le théorème de l'échantillonnage 
ou théorème de Shannon : 



Un signal qui ne contient pas de composantes à des fréquences supérieures ou 
égales à une valeur fm est entièrement déterminé par la suite de ses valeurs à des 



instants régulièrement espacés de la durée T = 



1 

2 fin ' 




0 fi . 1 f 

T f, = T 

Fig. 1 . 7 . Repliement de bande 



La fréquence d'échantillonnage d'un signal est ainsi déterminée par la limite 
supérieure de sa bande de fréquence. Dans la pratique on limite généralement par 
filtrage la bande du signal avant échantillonnage à la fréquence f e , à une valeur 
inférieure à /e/2 pour que le filtre de restitution soit réalisable. 

Il est intéressant de remarquer que la fréquence d'échantillonnage d'un signal 
est en fait déterminée par la largeur de bande qu'il occupe ; en effet la restitution a 
été illustrée sur la figure 1.6 pour un signal basse fréquence auquel a été associé un 
filtre passe-bas. On conçoit que le même raisonnement s’applique aussi à un signal 
occupant un domaine limité de l'espace des fréquences auquel serait associé un 
filtre passe-bande. Cette propriété est applicable en particulier aux signaux modu- 
lés et est utilisée dans certains types de filtres numériques. 

Le résultat donné à la fin du paragraphe 1.1.2 permet de présenter l’échan- 
tillonnage sous un autre aspect. En effet, la relation (1.57) montre que l’échan- 
tillonnage correspond à une décomposition du signal s(t) suivant l’ensemble des 



fonctions orthogonales 



sin n(t/T - n ) 
ic(r/T - n) 



et le théorème de Shannon exprime simple- 



ment la condition pour que cet ensemble forme une base de décomposition du 
signal. 



1.8 ÉCHANTILLONNAGE DE SIGNAUX SINUSOÏDAUX 
ET DE SIGNAUX ALÉATOIRES 

Les propriétés énoncées ci-dessus sont bien illustrées par l'échantillonnage de 
signaux sinusoïdaux, dont les particularités sont utilisables dans de nombreuses 
applications. 
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1.8.1 Signaux sinusoïdaux 

n 

Soit le signal s(t ) = cos (2k ft + cp), avec 0 =£ (p =s -, échantillonné avec la 
période T = l/f e = 1. 

Les échantillons sont donnés par la suite s (77) telle que : 

5 (n) = cos (2n fn + (p) 

Si le rapport f/f e = /est un nombre rationnel, il vient : 

/= N1/N2 avec NI et N2 entiers. 

Alors : 

s(n + N2) = cos [2tt f(n + N2) + cp] = s(n) 



La suite .s ( 77 ) présente la périodicité N2 et comprend au plus N2 nombres dif- 
férents. D'autre part la fréquence d'échantillonnage étant supérieure au double de 

1 

la fréquence du signal, on a nécessairement : N1/N2 < - . L'ensemble de N2 

échantillons différents permet de représenter un nombre de signaux sinusoïdaux 
égal au plus grand entier inférieur à N2/2. Par exemple si N2 = 8, avec l'ensemble 



des nombres : 2 cos \ 2n - + cpj, (n = 0, 1, ..., 7), il est possible de représenter les 



échantillons des 3 signaux sinusoïdaux : 



/ NI 

2 cos 1 2tt -g- t + 



avec NI = 1, 2, 3 




La figure 1.8 représente cet échantillonnage pour cp = 0; dans ce cas il suffit 
même de 4 nombres : ± 2 et ± V2. 
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Si l'on ajoute aux trois signaux sinusoïdaux de la figure 1.8, le signal continu 
de valeur 1 et le signal à la fréquence 1/2 d'amplitude 1 qui s’écrit : cos (nt), 
l'échantillonnage de cette somme donne des valeurs nulles, sauf aux instants mul- 
tiples de 8, où la valeur 8 est obtenue, comme le montre la figure 1.9. a. Le spectre 
de cette somme est obtenu directement en appliquant la relation : 

cos jc = - (e JX + e~’ x ). 

Il est formé de raies d'amplitude 1 aux fréquences multiples de 1/8 (fig. 1.9. b). Or 
ce spectre a été étudié au paragraphe 1.2, et l'on peut constater que la relation 
(1.27) se trouve vérifiée. 

La possibilité d'engendrer une gamme de signaux sinusoïdaux à partir d'un 
ensemble limité de nombres, stockés par exemple dans une mémoire, est utilisée 
dans les synthétiseurs de fréquence numériques. 




b) Spectre correspondant 



1 .8.2 Signaux aléatoires discrets 

Si le signal aléatoire s(t) est échantillonné avec la période supposée unitaire 
T = 1, il en résulte un signal aléatoire discret s(n), qui a par définition la même loi 
de probabilité de l'amplitude. Les résultats obtenus dans le cas continu se transpo- 
sent au cas discret, en particulier pour les signaux aléatoires stationnaires du 
second ordre ergodiques [5]. 
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Ainsi la fonction d'auto-corrélation du signal discret s(n ) est la suite r(n ) telle 
que : 

r(n) = E [.y (i ) . s (i - n)\ (1.58) 

C’est un échantillonnage de la fonction d'auto-corrélation r xx (x) du signal 
aléatoire continu définie par l'expression (1.34). Sa transformée de Fourier donne 
la densité spectrale énergétique O c (/) du signal discret, qui est liée à la densité 
spectrale O xx (/) du signal continu par une relation analogue à (1.55), c’est-à-dire : 

©e(/)= Y t) (L59) 

Si la fréquence d'échantillonnage n'a pas une valeur suffisante, ou si le spectre 
0 Vï (/) s’étend sur un domaine non borné, il y a repliement. 

L'hypothèse d’ergodicité pour le signal discret 5 (/z) conduit à la relation : 

r( - n) = 2NTÏ /= S n s(/) s(i ~' 7) (L60) 

Cette relation permet d'appliquer la notion de fonction d'auto-corrélation aux 
signaux déterministes. Ainsi, pour un signal périodique et de période N 0 , la fonc- 
tion d'auto-corrélation est la suite r(n) définie par : 

r(n) = Vf X s(i)s(i-n) (1.61) 

C’est une suite périodique et de même période. 

Exemple : 

s (n) = A sin ^2 ît j 



r{n) 



1 N o — t / i \ / i-n 

Nô ,?» A,Sin l 2 "N;J si "l 2 " Ni 



A 2 

r(n) = cos 



n 

2.71 XT 
N 0 



On retrouve bien la puissance du signal pour r (0) et la périodicité N 0 . 

Un signal aléatoire discret peut aussi être défini en tant que tel. Par exemple si 
la suite r(n) s'annule pour n + 0, le signal s(n) est un bruit blanc discret dont la 



densité spectrale est constante sur l'intervalle de fréquence 



1 

2 ’ 



1 

2 



. Ce signal a 



une réalité physique, c’est une suite de variables aléatoires non corrélées; pour 
l'obtenir il suffit de faire appel à un algorithme qui fournit des nombres statistique- 
ment indépendants. 
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1.8.3 Génération d'un bruit discret 

La génération de nombres aléatoires figure généralement au catalogue des fonc- 
tions des calculateurs scientifiques. Il est ainsi possible en logiciel de former une 
suite de nombres, utilisable comme signal de test en traitement numérique. 

Au paragraphe 1.8.1 on a montré qu’il est particulièrement simple de produire 
numériquement des signaux sinusoïdaux; de tels signaux peuvent aussi servir à 
simuler un bruit, par exemple par addition d’un grand nombre de sinusoïdes de fré- 
quences différentes, d’amplitude constante et de phase aléatoire ou pseudo-aléa- 
toire. Cette approche peut conduire à des réalisations particulièrement simples, 
comme la méthode qui a été normalisée pour l’appareillage de mesure utilisé dans 
les transmissions téléphoniques numériques. Cette méthode consiste à engendrer 
une séquence pseudo-aléatoire, qui est une suite périodique de 2 N - 1 bits compre- 
nant à une unité près autant de «zéros» que de «uns» et qui simule une suite aléa- 
toire dans laquelle les bits seraient indépendants et auraient la probabilité 1/2 de 
valoir «zéro» ou «un», ou pour centrer les variables, de valoir - 1/2 ou + 1/2. 

Si une opération de filtrage, qui en fait consiste en une sommation pondérée, 
est effectuée sur une telle suite, les nombres obtenus après filtrage suivent une loi 
de probabilité qui s’approche de la loi normale. 

Les séquences pseudo-aléatoires sont étudiées dans la référence [6], elles sont 
facilement obtenues à l’aide d’un registre à décalage à N bits, convenablement 
bouclé. La figure 1.10 donne un exemple, utilisé en appareillage de mesure, où 
N = 17. Le polynôme générateur s’écrit : 

g(jr) = 1 + x 3 + x 17 (1.62) 





Fig. 1.10. Générateur de séquence pseudo-aléatoire et loi de probabilité après filtrage 



la suite comprend 2 N - 1 = 131071 bits, elle est périodique et de période 
1 

T = (2 N - 1) .x, si t = -j- désigne la période de l’horloge du circuit. Le spectre est 
Jh 

formé de raies distantes de — . Pour / H = 370 kHz, l’espacement entre deux raies 
est de 2,8 Hz et l’on trouve 36 raies dans 100 Hz. 
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En opérant sur cette suite un filtrage à bande étroite qui ne conserve que la 
bande 450-550 Hz on obtient un signal approchant les caractéristiques gaussiennes, 
dont le facteur de crête est de 10,5 dB et qui constitue un excellent signal de test 
pour les équipements de transmission numérique. Si le filtrage est fait numérique- 
ment la suite de nombres obtenue peut être utilisée pour tester les équipements de 
traitement numérique. 



1.9 L'OPÉRATION DE QUANTIFICATION 

La quantification est l'approximation de chaque valeur du signal s(t) par un mul- 
tiple entier d'une quantité élémentaire q, appelée échelon de quantification. Si q 
est constant quelle que soit l'amplitude du signal, la quantification est dite uni- 
forme. Cette opération revient à faire passer le signal dans un organe qui possède 
une caractéristique en marche d'escalier, comme le montre la figure 1.11 pour 
q = 1, et fournit le signal sq(t). 




Fig. 1 . 11 . L’opération de quantification 

La manière dont l'approximation est faite définit le centrage de cette carac- 
téristique. Par exemple la figure représente le cas, appelé arrondi, où toute 
valeur du signal comprise entre ( n - l/2)q et ( n + H2)q est arrondie à nq. C’est 
l'approximation par défaut, qui est désignée, quand elle porte sur des nombres, 
par troncation et qui consiste à approcher par nq toute la valeur comprise entre 
nq et ( n + 1 )q ; la caractéristique se déplace alors de q!2 vers la droite sur l'axe 
des abscisses. 
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L'effet de cette approximation est de superposer au signal d'origine un signal 
d'erreur e(t) désigné par distorsion de quantification ou plus communément par 
bruit de quantification ; il vient : 

s(t)=s q (t) + e(t) (1.63) 

Une illustration est donnée par la figure 1.12 dans le cas de l'arrondi. Les 
amplitudes multiples impairs de q/2 sont appelées amplitudes de décision. 

L'amplitude du signal d'erreur est comprise entre - ql 2 et q/2. Sa puissance 
mesure la dégradation que subit le signal. 

Quand les variations du signal sont grandes par rapport à l'échelon de quanti- 
fication, c’est-à-dire que la quantification est faite avec suffisamment de finesse, le 
signal d'erreur est équivalent à un ensemble de signaux élémentaires, constitués 
chacun par un segment de droite (fig. 1.13). La puissance d'un tel signal élémen- 
taire de durée T s’écrit : 

i f 2 i / <7 \ 2 r 2 < j 2 

B = - \\t 2 dt = u (1.64) 

2 2 



S(t) 

S q (t) 




Fig. 1.12. Erreur de quantification 



Fig. 1.13. Signal d’erreur élémentaire 
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La valeur ainsi obtenue, B 



î! 

12 ’ 



est une estimation de la puissance du bruit de 



quantification suffisante dans la plupart des cas réels. 

La distribution spectrale du signal d'erreur est plus difficile à cerner. Le spectre 
du signal d'erreur élémentaire de la figure 1.13, E T (/), peut être calculé à partir de 
celui de sa dérivée. Ainsi en utilisant les relations (1.22) puis (1.12), on obtient : 



e,(/)=7^7?- 



sin {nfx) 
nfx 



cos {nfx) 



(1.65) 



Il apparaît que la plus grande partie de l'énergie se trouve au voisinage de la 
1 

fréquence Dans ces conditions la répartition spectrale du signal d'erreur dépend 

d'une part de la pente du signal élémentaire, c’est-à-dire en fait de la distribution 
statistique de la dérivée du signal s'{t), et d'autre part de la grandeur de l'échelon 
de quantification q par rapport au signal. La référence [7] donne le calcul de ce 
spectre pour un signal de bruit et fait apparaître un étalement en fonction de la fré- 
quence, quand le pas de quantification est suffisamment petit, sur un domaine qui 
couvre plusieurs centaines de fois la largeur de bande du signal. Si le signal à quan- 
tifier n’est pas un signal aléatoire, le spectre du signal d'erreur peut se concentrer 
sur certaines fréquences par exemple les harmoniques d'un signal sinusoïdal. 

Dans la conversion d'un signal analogique sous forme numérique, la quantifi- 
cation intervient conjointement avec l'échantillonnage, ces deux opérations étant 
réalisées successivement. Bien que l'échantillonnage soit en général fait en pre- 
mier, il est équivalent de faire la quantification d'abord et l'échantillonnage 
ensuite, à une fréquence f e habituellement un peu supérieure au double de la lar- 
geur de bande du signal. Dans ces conditions le signal d'erreur a souvent un 
spectre qui s’étend bien au-delà de la fréquence d'échantillonnage, et comme c’est 
en réalité la somme du signal et du signal d'erreur qui est échantillonnée, le phéno- 
mène de repliement du spectre intervient et la totalité de l’énergie du signal d'er- 
reur se retrouve dans la bande de fréquences [- f e l 2, f e l 2], La plupart du temps les 
conditions sont remplies pour que la densité spectrale énergétique du bruit de 
quantification soit constante et l'on retiendra le résultat suivant : 

Le bruit produit dans l’opération de quantification uniforme avec un échelon 
q a une puissance qui s’exprime en général par B = q 2 /12 et présente une réparti- 
tion spectrale constante dans la bande de fréquences [- f e l2,f e l2\. 

Il faut remarquer que la quantification des petits signaux, ceux dont l'ampli- 
tude est de l'ordre de grandeur de l'échelon q, dépend beaucoup du centrage de la 
caractéristique. Par exemple avec le centrage de la figure 1.11 un signal sinusoïdal 
d'amplitude inférieure à q!2 est totalement supprimé. Il est possible cependant de 
coder convenablement ces petits signaux en leur superposant un signal auxiliaire 
de grande amplitude qui est éliminé par la suite. 

Le codage d'un signal introduit ainsi une limitation pour les faibles amplitudes 
mais il impose également une borne aux fortes amplitudes. 
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Le signal échantillonné et quantifié en amplitude est représenté par une suite de 
nombres presque toujours sous forme binaire. Si chaque nombre compte N bits, 
le nombre maximum d'amplitudes quantifiées qu'il est possible de distinguer 
s’élève à 2 N. Alors la gamme des amplitudes qu'il est possible de coder est sou- 
mise à une double limitation : vers les faibles valeurs elle se trouve limitée par 
l'échelon de quantification q et vers les fortes valeurs par 2N. q. Toute ampli- 
tude qui dépasse cette valeur ne peut être représentée et il y a écrêtage du signal. 
Il s’en suit une dégradation, par exemple par distorsion harmonique si le signal 
est sinusoïdal. 

Si la gamme des amplitudes à coder couvre le domaine [- A m , + A m ], il vient : 

A„, = 2 N . q!2 (1.66) 

et d'autre part, avec l'arrondi, le signal d’erreur e(t) est tel que : 

|e(0|^A m .2-N 



On appelle puissance de crête d'un codeur la puissance du signal sinusoïdal 
ayant l'amplitude maximale admissible sans écrêtage, A m . Elle s’exprime par : 



P 



C 



1 

2 



" 2 N .q V 
2 



2 2N -\q 2 



La figure 1.14 représente ce signal avec le pas de quantification et les amplitudes 
de décision. 




On définit la dynamique du codage comme le rapport de cette puissance de 
crête à la puissance du bruit de quantification ; c’est en fait le rapport signal à bruit 
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maximal pour un signal sinusoïdal avec codage uniforme. Cette dynamique s’ex- 
prime par la formule suivante : 

P c /B = (S/B) max = 2 2N - 3 .12 = 3/2. 2 2N 

ou plus commodément en décibels : 

P f /B =6,02 N + 1,76 dB (1.67) 

Cette formule, d'une grande utilité pratique, relie le nombre de bits du codage 
à la plage des amplitudes qui peuvent être codées. 

Très souvent cependant le signal à coder n’est pas un signal sinusoïdal. Il est 
possible toutefois de se ramener à ce cas si, pour le signal à coder, est définie une 
puissance équivalente de crête, qui est alors prise comme puissance de crête du 
codeur. Le cas des signaux aléatoires gaussiens, est particulièrement important car 
ils représentent convenablement beaucoup de signaux rencontrés en pratique. Il 
faut alors positionner correctement l'amplitude maximale du codeur par rapport à 
l'amplitude du signal, de façon à ce que la distorsion introduite par écrêtage reste 
dans les limites imposées. 

En examinant le tableau donné en annexe 2, on peut remarquer que la proba- 
bilité pour que l'amplitude d'un signal de moyenne nulle et de puissance o 2 
dépasse 3,4 a est inférieure à ICC 3 . La figure 1.15 donne un exemple de codage 
avec o = q. Il apparaît que la probabilité d'écrêtage est inférieure à 5 . 1CL 4 avec les 
valeurs de paramètres choisies. 




Fig. 1 . 15 . Codage d’un signal gaussien 
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Finalement, pour obtenir le rapport signal à bruit maximal associé à un signal 
donné, il faut considérer le rapport signal à bruit de crête : 

(S/B) c = = 3 . 2 N 

et soustraire le facteur de crête F c . L'expression générale du rapport signal à bruit 
maximal pour un signal quelconque est donc la suivante : 

(S/B) max = 6,02 N + 4,77 - F,. dB (1.67bis) 

Ce résultat est utilisé non seulement dans la spécification des codeurs de 
signaux analogiques mais aussi dans le traitement numérique pour la détermina- 
tion des mémoires de données et le cadrage des nombres. 

La dynamique de codage, pour un nombre de bits donné, peut être considéra- 
blement augmentée si le codage est fait avec un échelon de quantification qui varie 
avec l'amplitude du signal; c’est le codage non linéaire. De nombreuses lois de 
variation peuvent être envisagées. Cependant il en est une qui est particulièrement 
importante puisqu'elle a été normalisée par l'Union Internationale des Télé- 
communications (UIT) pour le codage de signaux téléphonique dans les réseaux 
de télécommunications, c’est la loi segmentée à 13 segments [8]. 



1.11 CODAGE NON LINÉAIRE 

SUIVANT UNE LOI SEGMENTÉE 

Dans l'opération de codage non linéaire suivant la loi segmentée à 13 segments, les 
amplitudes positives et négatives à coder sont divisées en 7 plages, à chacune des- 
quelles est associé un échelon de quantification dont la grandeur résulte de la mul- 
tiplication d'un échelon élémentaire q par une puissance de 2. Cette opération 
peut être considérée comme résultant d'un codage linéaire précédé d'une com- 
pression selon laquelle le signal x est transformé en signal y conformément aux 
relations suivantes : 

1 + ln A \x | 1 

y = signe (x) . x + ln A P our A ^ M ^ 1 

, , ( 1 - 68 ) 

A \x\ 1 

y = signe (x ) . 1 + ln A pour 0 « |x| *£ — 

Le paramètre A détermine l'augmentation de la dynamique du codeur; la 
valeur retenue est A = 87,6. Finalement la caractéristique de compression suivant 
la loi A à 13 segments est donnée par la figure 1.16 et décrite comme suit : 

si 0 =£ \x\ =£ — , alors : y = 16x 
64 

1 . . 1 

64 1 1 32 



y = 8x + 1/8 
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1 

32 


1 

1 16 


y = 4x+ 1/4 


1 

Î6* 1 * 


1 

1 8 


y — 2x + 3/8 


1 

- ^ \x\ 
8 11 


« 1/4 


y — x + 1/2 


1 

- \x\ 

4 1 1 


« 1/2 


y = - je + 5/8 
2 


1 

2 11 


« 1 


V = \ x + 3/4 

4 




Fig. 1.16. Caractéristique de compression à 13 segments 

Cette caractéristique fait apparaître 7 segments de droite dans le quadrant 
positif et dans le quadrant négatif; les 2 segments qui entourent l'origine étant coli- 
néaires, au total la caractéristique compte bien 13 segments. 

La quantification des amplitudes y étant faite avec l'échelon q, celle des 
amplitudes x voisines de l'origine est faite avec l'échelon q/16, c’est-à-dire que la 
dynamique du codeur se trouve augmentée de 24 dB. Les amplitudes voisines de 
l'unité sont moins bien quantifiées puisque l'échelon se trouve multiplié par 4. La 
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puissance du bruit de quantification est fonction de l'amplitude du signal : pour 
chaque valeur il faut calculer un échelon moyen faisant intervenir la statistique du 
signal. 

La figure 1.17 donne le rapport signal à bruit en fonction du niveau du signal 
après codage, pour le codage à 8 bits linéaire et non linéaire d'un signal gaussien. 

Le niveau de référence pour le signal (0 dB) est la puissance de crête du 
codeur. On remarque l'extension de la dynamique due au codage non linéaire. 
Pour les amplitudes faibles la quantification correspond en fait à 12 bits. En réalité 
le signal codé suivant la loi non linéaire peut être obtenu à partir d'une quantifica- 
tion à 12 bits, suivie d'un traitement numérique qui est très proche de la conver- 
sion d'un nombre entier en un nombre à virgule flottante : 

Par exemple au nombre à 12 bits : +00010110110 

correspond le nombre à 8 bits : +100 0110 

par application de la loi de compression. 




Fig. 1.17. Codage à 8 bits linéaire et non linéaire d’un signal gaussien 

Les trois bits qui suivent le signe donnent le code du segment, ou exposant ; 
les quatre bits suivants indiquent la position sur le segment, ou mantisse. La dif- 
férence avec la conversion entier-virgule flottante apparaît au voisinage de l'ori- 
gine. 

Ce traitement nécessite pour sa mise en œuvre soit un réseau de portes c'est la 
réalisation parallèle, soit un registre à décalage associé à un compteur à 3 bits dans 
la réalisation série. On peut également utiliser une mémoire où est stockée la table 
de conversion. 

Une autre loi de codage non linéaire est également utilisée en télécommunica- 
tions, dite loi p à 15 segments. Elle correspond à la relation de compression sui- 
vante : 



y = signe (x) . 



In (1 + p|x|) 
ln (1 + p) 



pour - 1 x =£ 1 



(1.69) 



La valeur retenue pour le paramètre de compression est p = 255. 
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1.12 OPTIMISATION DU CODAGE 

En poursuivant dans la voie du perfectionnement du codage, si la densité de pro- 
babilité p (x) de l'amplitude du signal est connue, on peut déterminer la caractéris- 
tique de quantification qui, pour un nombre de bits N donné, minimise la puissance 
de la distorsion totale. 

Dans l'opération de quantification, la plage des amplitudes du signal est divisée 

M M 

en M = 2 N plages élémentaires x ; ) avec - — + 1 =£ i =£ — , et chaque plage 

élémentaire est représentée par une valeur y h comme le montre la figure 1.18. 
L'optimisation consiste à déterminer l'ensemble des valeurs x ; et y, qui minimise la 
puissance du signal d’erreur E 2 donnée par : 



M 




Fig. 1.18. Caractéristique de quantification optimale 



En prenant la dérivée par rapport aux variables x, et y t , on montre que l'opti- 
mum est obtenu si les relations suivantes sont vérifiées : 

1 MM 

*.■ = 2 (y< + + 1) p° ur - y + lss/ss y ~ 1 

[ x ‘ MM 

(x - y,) p (x) dx = 0 pour - — + 1 i =s — (1-70) 

1 L L 
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Ces relations permettent de déterminer la caractéristique de quantification. Si 
p (jc) est une fonction paire, on prend x 0 = 0 et l'on procède par itérations en choi- 

M 

sissant a priori une valeur de y v Si la relation (1.70) n'est pas satisfaite pour —, 
on reprend les calculs pour une autre valeur de y 1 et ainsi de suite [9]. 

Le tableau 1.1 donne la puissance du signal d'erreur E 2 obtenue avec un 
signal gaussien de puissance unité, pour différentes valeurs du nombre de bits 
N, d'une part avec le codage optimal, d'autre part avec un codage à échelon 
constant pour le meilleur cadrage de la caractéristique de quantification [9]. Le 
tableau 1.2, tiré de la référence [10], donne, dans les mêmes conditions, les 
valeurs correspondant à un signal à densité de probabilité exponentielle selon la 
relation (1.48). 



Tableau 1.1. - Codage d’un signal gaussien unitaire. 



N 

E, 2 


1 


2 


3 


4 


5 


Codage optimal 


0,3634 


0,1175 


0,03454 


0,0095 


0,0025 


Codage 
à échelon 
constant 


0,3634 


0,1188 


0,03744 


0,01154 


0,00349 


Entropie H 


1 


1,911 


2,825 


3,765 


4,730 



Tableau 1.2. - Codage d’un signal laplacien unitaire. 



N 

E 2 


1 


2 


3 


4 


5 


Codage optimal 


0,5 


0,1765 


0,0548 


0,0154 


0,00414 


Codage 
à échelon 
constant 


0,5 


0,1963 


0,0717 


0,0254 


0,0087 



L'optimisation du codage peut aussi être abordée sous l’aspect du contenu 
informationnel en introduisant la fonction entropie H définie par [2, Tome 3] : 

H = - X p, log 2 (p,) (1.71) 

i 



avec — — + 1 =£ i =£ — et où p, désigne la probabilité pour que l'amplitude du 
signal se trouve dans la plage représentée par la valeur y t . 
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Compte tenu de la relation : 

îi Pi = 1 

i 

l'entropie est nulle lorsque l'amplitude du signal se concentre sur une seule plage 
et elle est maximale lorsque cette amplitude est répartie uniformément; dans ce 
cas elle prend la valeur H max égale au nombre de bits N du codage : 

H max = log 2 (M) = N (1.72) 

où log 2 (M) représente le logarithme en base 2, ou binai, du nombre M. 

En fait l'entropie mesure l’écart d'une distribution de probabilité avec la dis- 
tribution uniforme. 

La caractéristique de quantification qui maximise l'entropie est donc celle qui 
conduit à des plages élémentaires correspondant à une distribution de probabilité 
uniforme. 

La dernière ligne du tableau 1.1 montre que pour un signal gaussien, le codage 
qui minimise la puissance du signal d'erreur conduit à des valeurs de l’entropie 
proches du maximum N. 



1.13 QUANTITÉ D'INFORMATION ET CAPACITÉ 
D'UN CANAL 



Les résultats obtenus sur l'échantillonnage et la quantification peuvent être utili- 
sés, à l'inverse, pour évaluer la quantité d'information portée par un signal ou pour 
déterminer la capacité d'un canal de transmission. 

Un canal réel de largeur de bande f m peut transporter 2 f m échantillons indé- 
pendants par seconde, comme le montre la figure 1.3, en remplaçant T par 2 f m . 
Quant à la quantité d'information par échantillon, elle dépend des puissances res- 
pectives du signal utile et du bruit, et de leurs distributions d'amplitude. Un cas 
particulier important est celui du canal gaussien [11]. 

Soit à transmettre un ensemble de M symboles de N bits chacun dans un canal 
en présence d'un bruit blanc gaussien de puissance o 2 h = B. 

Dans un hyperespace à M dimensions, les M symboles occupent le volume 
d'une hypersphère V M , défini par 



V M = 



'R 


c 


r M_ 1 d r 


... 


0 


; = i 



RM 

/(0,)d0, = — F(0) (1.73) 



En supposant une répartition uniforme des symboles dans l'hypersphère de 
rayon R, le signal correspondant a pour énergie : 



E ,= 



1 

v7 



MJ 



'R 


r 


r 2 r M-l d r 


... 


0 


i=i 



M-l 



m de, = R 2 



M 



M + 2 



(1.74) 
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La quantité d'information transmise est de MN bits. À chaque ensemble de 
bits possible, on peut associer dans l'hypersphère un volume V s égal à : 



V 



V = — — = — F(0)l — 

s OMN A K ' U ON 



1 



2 MN M 



R\. 



. M 



2 N j 



(1.75) 



À chaque ensemble de bits est associé un bruit à M composantes dont l'éner- 
gie s’écrit : E ; , = Mctg. Quand M tend vers l'infini, le point représentatif du bruit 
dans l'hypersphère se rapproche d'une sphère de rayon VMo 6 et centrée sur le point 
représentatif de l'ensemble des bits. En effet, pour M variables aléatoires gaussiennes, 



b(n ), la variable r 



M 



X b 2 {n) a pour moment d'ordre 1 : m 1 

n = 1 



M 2 
M+ 1 



et sa variance m 2 - m\ tend vers zéro quand M tend vers l'infini. Cette sphère a 
pour volume : 

(VMo fc ) M 



V fc = 



M 



-F(0) 



(1.76) 



Pour que la transmission se fasse sans erreur, il faut que cette sphère soit 
incluse dans le volume V s attribué à chaque ensemble de bits, c’est-à-dire : 



VMo ft < | 



(1.75) 



Or, quand M tend vers l'infini, d'après la relation (1.74), R 2 représente l'éner- 
gie de l'ensemble du signal, de puissance S et du bruit, soit 

R 2 = M(S + a g) (1.76) 

La condition de transmision sans erreur s’écrit alors 

S + al 

2 2N < — ^ (1.77) 

d'où : 

N<ilog2(l+^-) (1.78) 

Si le canal est réel et a une largeur de bande W et s'il est sans distorsion, les 
symboles peuvent être émis à la cadence 2W et la capacité asymptotique du canal 
s’écrit, en bits par seconde : 

C = Wlog 2(l+|) (1.79) 



Il faut bien noter qu'une telle capacité suppose : 

- un canal sans distorsion, 

- un bruit blanc gaussien, 

- un retard à la transmission infini. 
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En pratique, l'égalisation des canaux et le codage correcteur d'erreurs permet- 
tent de se rapprocher de cette limite avec un retard de transmission fini. 



1.14 LES REPRÉSENTATIONS BINAIRES 

Il existe diverses façons d'établir la correspondance entre l'ensemble des ampli- 
tudes quantifiées et l'ensemble des nombres binaires qui doivent les représenter. 
Les signaux à coder ayant des amplitudes en général positives et négatives, les 
représentations préférées sont celles qui conservent l'information de signe. Les 
plus courantes pour les codages à échelon constant sont les suivantes : 

- signe et valeur absolue 

- binaire décentré 

- complément à 1 

- complément à 2. 

Les définitions et particularités de ces représentations sont données dans la 
référence [12], le tableau 1.3 les définit pour 3 bits. 

Les représentations en signe et valeur absolue et en binaire décentré sont les 
plus commodes pour la conversion Analogique/Numérique; les deux autres sont 
surtout utilisées dans les circuits de calcul numérique; elles sont présentées en 
détail par la suite. 

Comme indiqué au paragraphe 1.11, qui décrit une représentation particulière 
importante du domaine des télécommunications, le codage non linéaire permet 
d'augmenter considérablement la dynamique. Les machines de traitement, et en 
particulier celles qui sont à usage général, utilisent souvent des représentations 
dites à virgule flottante où chaque nombre comporte trois parties : le bit de signe, 
la mantisse et l'exposant. La mantisse représente une partie fractionnaire et l'expo- 
sant la puissance d'un nombre de base ; par exemple, en base 10 : + 0,719 x 10 5 . 



Tableau 1.3. - Représentations binaires pour codage linéaire. 



Nombre 


Signe et valeur 
absolue 


Binaire 

décentré 


Complément 
à 1 


Complément 
à 2 


+ 3 


0 1 1 


111 


0 1 1 


011 


+ 2 


010 


110 


0 1 0 


010 


+ 1 


00 1 


101 


00 1 


00 1 


+ 0 


000 


100 


000 


000 


-0 


100 


- 


1 1 1 


- 


-1 


101 


011 


1 1 0 


1 1 1 


-2 


110 


010 


101 


110 


-3 


1 1 1 


001 


100 


101 






(0 0 0) 




(10 0) 
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L'extension de la dynamique provient de l'effet multiplicatif introduit par l’ex- 
posant. Ainsi en base 2, pour un nombre à 6 bits d'exposant et 16 bits de mantisse, 
la dynamique correspond à 2 64 x 2 16 = 2 80 = 10 24 , soit 24 chiffres décimaux. Un 
grain supplémentaire est obtenu en prenant une base qui est elle même une puis- 
sance de deux, comme 8 ou 16, correspondant aux numérations octale ou hexadé- 
cimale respectivement. 

La présentation à virgule flottante entraîne cependant une complication des 
opérations arithmétiques et des circuits. 

Les nombres issus du codage se présentent, suivant la technique utilisée, soit 
sous forme parallèle, c’est-à-dire que les N bits sont disponibles sur N points de 
connexion en même temps, soit sous forme série, c’est-à-dire que les N bits appa- 
raissent successivement sur le même point de connexion, le signe d'abord et 
ensuite les bits de poids décroissants. La référence [13] décrit les principales tech- 
niques de conversion Analogique/Numérique. 



ANNEXE 1 : La fonction l(%) 



Suite des valeurs : I (n) = 



n 

sm|7t — 

20 , 



n 

20 



pour 0 « n « 159 avec n = k + 20 N. 



k 


N = 0 


N = 1 


N = 2 


N = 3 


N = 4 


N = 5 


N = 6 


N = 7 


0 


1 


0 


0 


0 


0 


0 


0 


0 


1 


0,99589 


- 0,04742 


0,02429 


- 0,01633 


0,01229 


- 0,00986 


0,00823 


- 0,00706 


2 


0,98363 


- 0,08942 


0,04684 


- 0,03173 


0,02399 


- 0,01929 


0,01613 


- 0,01385 


3 


0,96340 


- 0,12566 


0,06721 


- 0,04588 


0,03482 


- 0,02806 


0,02350 


- 0,02021 


4 


0,93549 


- 0,15591 


0,08504 


- 0,05847 


0,04455 


- 0,03598 


0,03018 


- 0,02599 


5 


0,90032 


- 0,18006 


0,10004 


- 0,06926 


0,05296 


- 0,04287 


0,03601 


- 0,03105 


6 


0,85839 


- 0,19809 


0,11196 


- 0,07804 


0,05989 


- 0,04850 


0,04088 


- 0,03528 


7 


0,81033 


- 0,21009 


0,12069 


- 0,08466 


0,06520 


- 0,05301 


0,04466 


- 0,03859 


8 


0,75683 


- 0,21624 


0,12614 


- 0,08904 


0,06880 


- 0,05606 


0,04730 


- 0,04091 


9 


0,69865 


- 0,21682 


0,12832 


- 0,09113 


0,07065 


- 0,05769 


0,04874 


- 0,04220 


10 


0,63662 


- 0,21221 


0,12732 


- 0,09095 


0,07074 


- 0,05787 


0,04897 


- 0,04244 


11 


0,57162 


- 0,20283 


0,12329 


- 0,08856 


0,06910 


- 0,05665 


0,04800 


- 0,04164 


12 


0,50455 


- 0,18921 


0,11643 


- 0,08409 


0,06581 


- 0,05406 


0,04587 


- 0,03983 


13 


0,43633 


- 0,17189 


0,10702 


- 0,07770 


0,06099 


- 0,05020 


0,04265 


- 0,03707 


14 


0,36788 


- 0,15148 


0,09538 


- 0,06960 


0,05479 


- 0,04518 


0,03844 


- 0,03344 


15 


0,30011 


- 0,12862 


0,08185 


- 0,06002 


0,04739 


- 0,03914 


0,03335 


- 0,02904 


16 


0,23387 


- 0,10394 


0,06682 


- 0,04924 


0,03898 


- 0,03298 


0,02751 


- 0,02399 


17 


0,17001 


- 0,07811 


0,05071 


- 0,03753 


0,02980 


- 0,02470 


0,02110 


- 0,01841 


18 


0,10929 


- 0,05177 


0,03392 


- 0,02522 


0,02007 


- 0,01667 


0,01426 


- 0,01245 


19 


0,05242 


- 0,02554 


0,01688 


- 0,01261 


0,01006 


- 0,00837 


0,00716 


- 0,00626 
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ANNEXE 2 : La loi Normale Réduite 



X 2 


il 


* 2 

e 2 dx 

i 


2 

VXt . 


* 2 

e 2 dx 

i 


X 


10 5 -/ (x) 


100 P 


X 


X 


100 P 


0 


39894 


100 


0 


0 


100 


0,2 


39104 


95 


0,0627 


0,2 


84,148 


0,4 


36827 


90 


0,1257 


0,4 


68,916 


0,6 


33322 


85 


0,1891 


0,6 


54,851 


0,8 


28969 


80 


0,2533 


0,8 


42,371 


1 


24197 


75 


0,3186 


1 


31,731 


1,2 


19419 


70 


0,3853 


1,2 


23,014 


1,4 


14973 


65 


0,4538 


1,4 


16,151 


1,6 


11092 


60 


0,5244 


1,6 


10,960 


1,8 


7895 


55 


0,5978 


1,8 


7,186 


2 


5399 


50 


0,6745 


2 


4,550 


2,2 


3547 


45 


0,7554 


2,2 


2,781 


2,4 


2239 


40 


0,8416 


2,4 


1,640 


2,6 


1358 


35 


0,9346 


2,6 


0,932 


2,8 


792 


30 


1,0364 


2,8 


0,511 


3 


443 


25 


1,1503 


3 


0,270 


3,2 


238 


20 


1,2816 


3,2 


0,137 


3,4 


123 


15 


1,4395 


3,4 


0,067 


3,6 


61 


10 


1,6449 


3,6 


0,032 


3,8 


29 


5 


1,9600 


3,8 


0,014 


4 


13 


1 


2,5758 


4 


0,006 


4,2 


5,9 


0,1 


3,2905 


4,5 


0,00068 


4,4 


2,5 


0,01 


3,8906 


5 


0,000057 


4,6 


1 


0,001 


4,4172 


5,5 


0,000004 






0,0001 


4,8916 










0,00001 


5,3267 







Approximation pour les grandes valeurs du paramètre À. : 



P = 



3 

J 



1 

X 



e 



x. 

T 
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20 log A 




Fig. 1.19. Loi Normale Réduite. Courbe donnant P en fonction de 20 log X 
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EXERCICES 

1 Soit le développement en série de Fourier de la Fonction i{t) périodique et de 

X X 

période T, nulle sur toute la période, sauf dans l’intervalle - - =£ t =s - où elle prend la 
la valeur 1. 

T T 

Donner la valeur des coefficients pour x = - et x = - . 

F 2 3 

Vérifier que le développement conduit bien à : 7(0) = 1 et tracer la fonction quand le 
développement est limité à 5 termes. 

2 Analyser l'échantillonnage à la fréquence fe du signal 5 (t) = sin (nf e t + (p) quand (p 



varie de 0 à - . 

2 

Examiner la reconstitution de ce signal à partir des échantillons. 

3 Calculer la distorsion d’amplitude apportée à un signal restitué par des impulsions 
dont la largeur est égale à la moitié de la période d’échantillonnage. 

4 Echantillonnage passe-bande d’un signal occupant la bande de fréquence [ f - f 2 ]. 
Quelles sont les conditions à imposer à la fréquence f pour que ce signal puisse être échan- 
tillonné directement à une fréquence comprise entre f 2 et 2 f 2 . 

5 Analyser l’échantillonnage du signal suivant : 

n t 



,,(0= Z s in(2tt 8T 



et le comparer à celui du signal : 

Int 

s r (t)=Xœ S (2n -- 

Montrer par une étude des spectres que l’ensemble des deux suites d’échantillons constitue 
l’échantillonnage d’un signal complexe. 

6 Soit s ( t ) le signal défini par l’égalité : 
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Ce signal est échantillonné avec la période T = 1. Quelle est la valeur maximale que peut 
prendre s(n) avec n entier? Montrer qu’il existe un ensemble de valeurs <p k ( k = 1, 2, 3, 4) 
qui minimisent la valeur maximale des s (n). Peut-on généraliser cette propriété ? 

7 Un synthétiseur de fréquence numérique est construit à partir d’une mémoire morte 
de 16A:bits ayant un temps d’accès de 500 ns. Sachant que les nombres qui représentent les 
échantillons de signaux sinusoïdaux comptent 8 bits, quelles sont les caractéristiques du syn- 
thétiseur, gamme et pas de fréquence, qui peuvent être obtenues ? 

8 Quelle est la loi de probabilité des amplitudes du signal sinusoïdal suivant : 

/ t 

s(t) = A cos I 2 ji - 

Calculer la fonction d’autocorrélation correspondante. 

Donner la fonction d’autocorrélation d’un signal aléatoire gaussien stationnaire dont le 
spectre a une répartition uniforme dans la bande de fréquence [/i,/ 2 ]. 

9 Calculer le spectre d’une suite d’impulsions de largeur T/2, séparées de T, chaque 
impulsion ayant la probabilité p d’apparaître. Examiner en particulier le cas où p = 1/2. 

Que devient ce spectre si ces impulsions constituent une séquence pseudo-aléatoire de 
longueur 2 4 - 1 = 15 engendrée par un registre à 4 bits, suivant le polynôme g(x) = x 4 + x + 1 ? 

10 Un signal sinusoïdal à la fréquence 1 050 Hz est échantillonné à 8 kHz et codé à 
10 bits. Quelle est la valeur du rapport signal à bruit maximal? Quel est le niveau par rap- 
port au signal du bruit de quantification mesuré dans la bande de fréquence 300-500 Hz ? 
Même question si la fréquence d’échantillonnage est portée à 16 kHz. 

t \ Tt 

2n; - + cp avec 0 =s (p - est échantillonné avec la période T = 1 et 
8 / 2 

codé à 5 bits. 

Dans le cas où cp = 0 calculer la puissance et le spectre du bruit de quantification. 

Comment évolue ce spectre en fonction de la phase cp? 

12 Soit une échelle de codage où l’échelon a la valeur q. Étudier la quantification du 
signal i 1 (t) = a . ç.sin (<n h t) pour - 1 =s a =s 1, en fonction du centrage de la caractéristique 
de quantification. Donner l’enveloppe du signal restitué après décodage et filtrage étroit 
autour de la fréquence ci),. 

Au signal Sj (f) on superpose le signal s 2 (f) = 10 . q . sin c o 2 t. Étudier l’enveloppe du 
signal restitué dans ces conditions. 

13 Soit à coder un signal gaussien. Combien de bits faut-il pour que le rapport signal à 
bruit de quantification soit meilleur que 50 dB ? Peut-on réduire ce nombre si l’on admet un 
écrêtage pendant 1 % du temps. 

14 Le signal s(t) = A sin (2n . 810.f) est codé à 8 bits. L’échelon de quantification 
ayant la valeur q , tracer la courbe donnant le rapport signal à bruit de quantification en 
fonction de l’amplitude A lorsque cette amplitude varie de q al 1 .q. Même question lorque 
le codage est du type non linéaire suivant la loi à 13 segments. 

15 Calculer les limites des plages d’amplitude élémentaires pour le codage optimal à 
2 bits d'un signal gaussien de puissance unité. 



11 Le signal sin 




Chapitre 2 



La transformation de Fourier 



La transformation de Fourier Discrète s'introduit quand il s'agit de calculer la 
transformée de Fourier d'une fonction à l'aide d'un calculateur numérique. En 
effet un tel opérateur ne peut traiter que des nombres et de plus en quantité limi- 
tée par la taille de sa mémoire. Il s’en suit que la transformée de Fourier : 



S (/) = 



s (t) e~> 2 n f‘ dt 



doit être adaptée, d'une part en remplaçant le signal s(t ) par des nombres s (V/T) 
qui représentent un échantillonnage de ce signal et d'autre part en limitant l’en- 
semble des nombres sur lesquels portent les calculs à une valeur finie N. Le calcul 
fournit alors des nombres S* (/) définis par : 



S *(/)= X s(nT) 

n = 0 



e -j2nfnT 



Comme le calculateur est limité dans sa puissance de calcul, il ne peut fournir ces 
résultats que pour un nombre limité de valeurs de la fréquence /, qu'il est naturel 
de choisir multiples d'un certain pas de fréquence A f Alors : 

N - 1 

S*(kAf)= X s(nT)ed 2nnkA f 1 ' 

n = 0 



Les conditions dans lesquelles les valeurs calculées constituent une bonne approxi- 
mation des valeurs recherchées sont étudiées par la suite. Un choix simplificateur 

intéressant consiste à prendre : A / = . On peut alors vérifier qu'il existe seule- 

ment N valeurs différentes dans la suite des S* (/(/NT), qui est une suite périodique 
et de période N puisque 



S* [(k + N)/NT] = S* (U/NT ) 
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D'autre part la transformée ainsi calculée se présente sous la forme de valeurs dis- 
crètes et d'après le paragraphe 1.6 sur l'échantillonnage en fréquence, cette pro- 
priété est caractéristique du spectre des fonctions périodiques. On peut donc consi- 
dérer que la suite des S* (/c/NT ) est obtenue par transformation de Fourier de la 
suite des s (nT ) qui est une suite périodique et de période NT. 

La transformation de Fourier discrète (T.F.D.) et la transformée inverse éta- 
blissent des relations entre ces deux suites périodiques. 

La définition, les propriétés, les techniques de calcul et les applications de la 
T.F.D. ont fait l'objet de nombreux articles et ouvrages, parmi lesquels on peut 
citer les références [1, 2, 3, 4], 



2.1 DÉFINITION ET PROPRIÉTÉS DE LA TFD 



Soit deux suites de nombres complexes x(n) et X(/c), périodiques et de période N. 
La transformée de Fourier Discrète et la transformée inverse établissent entre ces 
deux suites les relations suivantes respectivement : 



-| N - 1 nk 

X{k) = — Z x(n) e l2K ~H 
N n = 0 


(2.1) 


N - 1 kn 

x(n) = Z X(k)e’ 2n n 

* = 0 


(2.2) 



La position du facteur d’échelle 1/N est choisie pour que les X(/c) soient les 
coefficients du développement en série de Fourier de la suite x (n). Cette transfor- 
mation possède les propriétés suivantes : 

- Linéarité : si x ( n ) et y (n) sont deux suites de même période, dont les transfor- 
mées sont X(/c) et Y (k) respectivement, la suite de même période v(n) = x(n) + 
Ky (n) où X est un scalaire a pour transformée : 

V(/c) =X(k) + X.Y{k) 

- Une translation des x(n) entraîne une rotation de phase des X (k) : En effet 
calculons la transformée X n .(k) de la suite x(n - n 0 ). 

N - 1 nk 

X (k) = Z x(n-n 0 )e n = X(k)e^i 2ïïn o klN 

0 n = 0 



Une translation de x(n) de n 0 entraîne sur X(/c) une rotation de la phase d'un angle 

. , , 0 «o k 
égal a : 2n 



- Symétrie : si la suite x (n) est réelle les nombres X(/c) et X(N - k) sont com- 
plexes conjugués : 

N n (N - k) 

X(N -k)= Z x(n)e' 2K ~^~ =X(/c) 

;i = 0 




